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ABSTRACT 

This  investigation  demonstrates  how  the  maximum  moisture  absorbed 
(that  is  the  "wet”  condition)  in  [0/±45/90]j  and  [0/90/-453g  laminates 
fabricated  from  T300/520a  significantly  alters  the  dry  stress  state  and 
subsequent  damage  development  along  the  laminate  free  edge. 

Emphasis  is  placed  on  using  reasonable  approximations  for  wet,  dry, 
and  out-of-plane  elastic  properties  since  these  properties  are 

required  to  predict  the  damage  free  stress  state  at  the  laminate  edge. 
Classical  laminate  theory  and  a  finite  element  model  were  used  to 
predict  stress  states  prior  to  the  first  formation  of  damage.  Crack 
patterns  characteristic  of  the  laminate  in  a  wet  or  dry  condition  were 
also  predicted  using  a  shear  lag  model.  Development  of  edge  damage  was 
recorded  and  observed  during  the  test  by  transferring  an  image  of  the 
damage  from  the  edge  surface  on  to  a  thin  acetate  sheet  (replica 
technique),  such  that  the  damage  imprinted  on  the  acetate  sheet  could 
be  immediately  viewed  on  a  microfiche  card  reader. 

Moisture  was  shown  to  significantly  alter  the  interior  and  edge 
dry  stress  states  due  to  swelling  and  a  reduction  of  elastic 
properties.  Moisture  also  reduces  the  transverse  strength  in  the 
90°  olies  such  that  the  first  formation  of  damage  in  a  wet  [0/±45/90]^ 


laminate  is  a  simultaneous  occurrence  of  delaminations  and  transverse 
cracks  in  the  90°  plies.  A  model  was  developed  in  order  to  predict 
changes  in  first  ply  failure  laminate  loads  due  to  differences  in 
stacking  sequence  together  with  a  wet  or  dry  environmental  condition. 

Although  moisture  was  shown  to  significantly  alter  the  first 
formation  of  damage,  the  crack  patterns  prior  to  fracture  were  not 
significantly  altered  by  moisture  absorption.  Consequently,  differences 
between  wet  and  dry  laminate  static  or  residual  experimental  strengths 


were  small. 
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I.  INTRODUCTION 


Recent  emphasis  on  environmental  degradation  of  graphite/epoxy 
fiber- reinforced  composites  is  due  to  the  increased  use  of  these  high 
performance  materials  in  aerospace  structural  applications.  It  is 
well  established  that  atmospheric  moisture  which  is  absorbed  by 
diffusion  into  the  epoxy  matrix  degrades  those  lamina  properties 
which  are  matrix  dependent  [1].  The  accumulative  moisture  absorbed 
by  T300/5208  graphite/epoxy  when  exposed  to  in-service  environments 
[2]  is  significant  in  terms  of  reduced  lamina  properties. 

Hygrothermal  degradation  of  graphite/epoxy  could  be  attributed  to 
degradation  of  the  fiber,  matrix,  or  fiber/matrix  interface.  It  is 
generally  accepted  that  fiber  properties  are  unaffected  by  moisture 
[1]  since  moisture  has  little  effect  on  lamina  properties  which  are 
fiber  dominated.  The  fiber/matrix  interface  strength  is  reduced  due 
to  moisture  [3].  This  can  be  attributed  to  a  combination  of  fiber/ 
matrix  chemistry  [4]  (i.e.  graphite  fiber  sizing)  and  a  residual  stress 
state  at  the  graphite/epoxy  interface  [5].  It  is  therefore  generally 
accepted  that  most  cracks  which  result  from  a  mechanically  applied 
load  in  graphite/epoxy  materials  initiate  at  the  fiber/matrix  interface 
When  moisture  is  absorbed  overall  degradation  of  the  matrix  exists 
since  most  epoxy  resins  are  susceptible  to  plasticization, 
enhanced  viscoelastic  response,  together  with  a  reduction  in  glass 
transition  temperature,  ultimate  strength,  and  stiffness  properties 
[1,6]. 

Degradation  of  the  epoxy  matrix  and  fiber/matrix  interface  is 
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the  result  of  volumetric  diffusion  [7]  of  water  molecules  which  attach 
themselves  as  hydrogen  bonded  molecules  onto  the  long  epoxy  polymer 
chains.  This  diffusion  increases  the  epoxy  "free  volume"  [1]  which 
results  in  swelling.  The  rate  of  moisture  absorbtion  can  be  conven¬ 
iently  accelerated  by  exposure  at  elevated  temperatures.  Unfortunately 
accelerated  moisture  absorption  will  produce  matrix  cracking  if  the 
temperature  exceeds  the  glass  transition  temperature  [8].  This  matrix 
damage  is  usually  near  the  surface  and  is  attributed  to  a  combination 
of  matrix  plasticization  and  a  residual  stress  state  which  is  created 
when  a  large  gradient  in  rxiisture  concentration  profiles  causes  the 
dry  surface  to  shrink  upon  desorption  [9].  Although  there  are  many 
more  interesting  damage  mechanisms  which  explain  the  formation  of 
cracks  in  epoxy  resins,  many  of  these  mechanisms  are  worst  cases  of 
laboratory  induced  degradations.  Once  these  worst  case  damage 
mechanisms  are  understood  they  are  usually  eliminated  from  a  materials 
application  viewpoint. 

Although  it  is  instructive  to  study  these  worst  case  mechanis.as , 
the  emphasis  in  this  study  is  to  choose  a  material  system  which 
minimizes  these  worst  case  cracking  events.  Eliminating  these  worst 
case  events  results  in  a  less  com.plex  model  which  can  then  be  used  to 
explain  how  damage  developes  in  an  environmentally  conditioned  lamin¬ 
ate  when  a  mechanical  load  is  applied.  Crossman  [10]  demonstrated 
that  the  strength  and  elastic  properties  of  quasi -isotropic  laminates 
fabricated  from  T300/5208  were  reduced  when  moisture  was  absorbed  (that 
is,  the  "wet"  condition),  and  that  no  damage  resulted  from  the 
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dbsorption  of  moisture.  Crossmdn  also  showed  that  the  viscoelastic 
response  of  the  wet  T300/52C8  laminates  was  negligible  when  com¬ 
pared  to  the  viscoelastic  response  of  wet  T300/5209  laminates.  This 
is  partially  due  to  the  lower  glass  transition  temperature  of  T300/5209. 
In  summary  quasi-isotropic  laminates  fabricated  from  T30C/5208  will 
behave  elastical ly  in  either  the  wet  or  dry  scate  such  that  wet  lamina 
strength  and  elastic  properties  are  lower  than  the  dry  properties  and 
no  cracks  result  when  T300/5208  laminates  absorb  moisture.  It  follows 
that  the  formation  of  any  cracks  in  a  quasi-isotropic  laminate  fab¬ 
ricated  from  T300/5208  will  be  the  result  of  an  applied  mechanical 
load  acting  together  with  the  residual  hygro-thermal  stress  state. 

Since  wet  or  dry  quasi-isotropic  laminates  fabricated  from  T300/ 
5208  behave  elastically,  the  interior  in-plane  stress  state  existing 
prior  to  the  initial  formation  of  damage  can  be  calculated  from 
classical  laminate  plate  analysis  using  wet  or  dry  elastic  properties. 
Using  a  laminated  plate  analysis  Kim  and  Hahn  [11]  predicted  the  stress 
state  in  a  wet  and  dry  [0/+45/-45/90]^  laminate  fabricated  from  T300/ 
5208  prior  to  formation  of  the  first  90°  ply  crack.  Good  correlation 
between  experimental  and  predicted  first  90°  ply  failure  stress  in 
both  wet  and  dry  conditions  was  demonstrated  by  using  averaged  wet  and 
dry  elastic  lamina  properties.  In  Kim's  model  the  wet  residual  stress 
state  was  assumed  zero  after  absorbing  1.3  percent  moisture  (by  weight 
gained)  which  implies  that  all  damage  is  due  to  a  state  of  stress 
resulting  only  from  an  applied  load.  The  first  ply  failure  loads 
were  predicted  using  a  stress  failure  criterion  along  with  the  laminate 


stress  state  predicted  by  laminated  plate  analysis  using  constant 
coefficients  of  expansion. 

The  stress  state  in  finite  width  quasi-isotropic  laminates  is  not 
uniform  through  the  width  as  assumed  by  classical  laminated  plate 
analysis.  Three  dimensional  stress  analysis  demonstrates  that  the 
stress  components  which  are  acting  perpendicular  to  the  thin  quasi¬ 
isotropic  laminate  plane  are  negligible  away  from  the  free  edge  but 
can  exceed  ply  strengths  within  a  thin  boundary  layer  near  the  free 
edge  [10].  These  out-of-plane  stresses  are  the  result  of  a  mismatch 
of  Poisson's  ratios  and  coefficients  of  thermal  expansion  for  each 
layer  in  the  laminate  when  loaded  mechanically  or  thermally.  The 
magnitude  of  the  out-of-plane  tensile  stresses  near  the  free  edge 
can  cause  interply  cracks  [12]  which  are  called  delamiriations .  No 
delaminations  occur  when  the  stacking  sequence  is  altered  to  give 
compressive  out-of-plane  stresses.  Depending  on  the  stacking  sequence 
the  damage  which  develops  along  the  free  edge  of  the  laminate  when 
loaded  can  be  a  combination  of  delaminations  between  layers  or 
transverse  cracks  within  layers.  Previous  studies  [13]  have  demon¬ 
strated  that  the  damage  which  develops  along  the  free  edge  of  a 
[0/+45/-45/90]^  laminate  (type  I)  is  entirely  different  from  the  free 
edge  damage  for  a  [0/90/+45/-45]g  laminate  (type  II).  Differences  in 
damage  states  observed  along  the  free  edge  can  ultimately  influence  the 
final  laminate  strength.  For  type  I  and  type  II  laminates  fabricated 
from  T300/5208  there  can  be  as  much  as  30  percent  difference  in  dry 
laminate  strengths  [14,15].  Tlierefore,  when  considering  damage  leading 


to  final  lamin.'.te  failure,  the  emphasis  changes  from  individual  trans¬ 
verse  cracks  to  the  development  of  an  entire  damage  state  near  the 
free  edge  prior  to  failure.  In  this  investigation  only  damage  at 
the  laminate  free  edge  is  investigated.  No  attempt  was  made  to 

investigate  damage  away  from  the  free  edge- 

Due  to  the  heterogeneous  nature  of  fiber-reinforced  laminates  the 
damage  which  develops  in  a  laminate  when  loaded  is  a  composition  of 
interply  delamination  cracks  and  cracks  within  individual  plies  which 
grow  transverse,  longitudinal,  or  at  an  angle  to  the  load  axis.  Unlike 
the  single  critical  crack  in  homogeneous  materials,  each  crack  which 
exists  in  a  heterogeneous  laminate  cannot  be  evaluated  as  an  isolated 
event  which  grows  in  a  self-similar  fashion.  Instead,  each  crack  in 
a  heterogeneous  laminate  is  a  component  of  a  damage  state,  and  the 
laminate  response  is  influenced  more  by  the  development  of  this  damage 

state  than  by  the  behavior  of  a  single  crack. 

Based  on  extensive  experimental  data,  Reifsnider  et  al  [16]  have 
shown  that  different  crack  types  consistently  develop  into  a  stable 
pattern  o-  "damage  state"  which  is  characteristic  of  the  laminate.  A 
rational  mechanistic  approach  can  be  used  to  predict  a  stable  pattern 
of  cracks  in  the  off-axis  plies  prior  to  laminate  failure  [17],  This 
characteristic  damage  state  (CDS)  could  be  used  to  define  the  stress 
and  state  of  strength  prior  to  laminate  failure.  Recent  observations 
[17]  strongly  suggest  that  the  CDS  could  be  independent  of  load 
history  and  that  the  CDS  is  a  laminate  property.  In  summary,  the 
philosophy  of  CDS  is  best  stated  by  Reifsnider  and  Masters  [17], 
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"From  the  standpoint  of  mechanics,  the  COS  has  the  same  significance 
as  the  single  crack  for  homogeneous  materials  ir:  the  sense  that  it 
is  the  well-defined  damaged  physical  state  from  which  the  fracture 
event  develops." 

As  already  discussed,  previous  investigations  have  shown  that  the 
formation  of  damage  in  quasi-isotropic  laminates  depends  on  material, 
stacking  sequence,  residual  cu**ing  stresses,  and  environmental  condi¬ 
tioning.  In  particular  the  present  investigation  is  primarily 
concerned  with  the  effect  of  residual  cure  stresses,  swelling  due  to 
moisture  absorption,  and  mechanical  loads  on  tne  CDS  in  type  I  and 
type  II  laminates  fabricated  from  T300/5208  graphite/epoxy.  The 
objectives  of  this  study  are  to  initially  isolate  unique  free  edge 
damage  states  in  type  I  and  type  II  laminates  and  show  how  these  unique 
damage  states  develop  into  the  CDS  when  the  laminate  is  mechanically 
loaded. 

Unique  free-edge  damage  states  in  quasi-isotropic  laminates  are 
obtained  experimentally  by  environmentally  conditioning  type  I  and 
type  II  laminates.  The  differences  in  the  laminate  wet  and  dry  damage 
free  stress  states  uniquely  influence  the  free  edge  damage  which 
develops  when  the  laminates  are  mechanically  loaded.  While  the  laminate 
load  is  held  constant,  damage  along  the  free  edge  is  recorded  by 
replicating  an  image  of  damage  from  the  free  edge  surface  on  to  the 
surface  of  an  acetate  strip.  Analytic  models  wnich  predict  the  state 
of  stress  prior  to  the  formation  of  damage  and  the  characteristic 
spacing  of  ply  cracks  along  the  free  edge  are  compared  with  the  damage 
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recorded  on  the  replicas  while  the  laminates  v^ere  cyclically  loaded 
or  incrementally  loaded  to  failure. 

If  the  initial  damage  free  stress  state  existing  near  the  free 
edge  uniquely  influences  the  subsequent  damage  state  and  laminate 
strength  then  accurate  out-of-plane  elastic  properties  must  be  used 
when  evaluating  this  stress  state.  To  date  only  Dean  and  Turner  [18], 
Ishikawa,  Koyama,  and  Kobayashi  [19],  and  Kriz  and  Stinchcomb  [20]  have 
obtained  reasonable  estimates  of  out-of-plane  lamina  properties.  As 
pointed  out  by  Crossman  [14],  out-of-plane  properties  used  in  most 
free  edge  stress  analyses  are  rough  approximatiors;  therefore,  only 
trends  in  stress  fields  can  be  demonstrated.  In  this  investigation 
accurate  wet,  and  dry  out-of-plane  elastic  lamina  properties  will  be 
used  in  various  stress  analysis  methods  to  evaluate  the  wet  and  dry 
stress  state  both  near  and  av/ay  from  the  free  edge. 


II.  EXPERIMENTAL  PROCEDURE 


2.1  Environmental  Conditioning 

The  specimens  listed  in  Table  1  were  fabricated  by  Southwest  Research 
Institute  and  McDonnell  Douglas  using  the  same  reconmended  procedure 
for  curing  Narmco  T300/5208.  Half  of  the  specimens  were  exposed  to  05?i 
RH  at  70°C,  and  the  remaining  specimens  were  dried  in  a  dry  nitrogen 
gas  oven  at  GB^C.  Although  vacuum  ovens  have  been  commonly  used  for 
moisture  desorption  at  elevated  temperatures,  damage  is  minimized 
when  dry  nitrogen  gas  is  used  [21],  Moisture  absorbed  or  lost  was 
measured  as  a  change  in  weight.  All  specimens  labeled  as, "WET"  were 
in  a  condition  of  maximum  absorbed  moisture  which  occurr'ed  when  no 
additional  increase  in  specimen  weight  could  be  measured.  All  specimens 
labeled  as  "DRY"  were  in  a  condition  of  total  absence  of  diffused 
moisture  which  was  obtained  when  no  additional  loss  in  specimen  weight 
could  be  measured. 

2.2  Unidirectional  Tension  Tests 

Both  wet  and  dry  [Og]  and  [90g]  specimens  listed  in  Table  1  were 
loaded  to  failure  in  tension  using  an  Instron  load  frame  with  a 
crosshead  speed  of  0.1  inches  per  minute.  Biaxial  strains  were  measured 
at  the  center  of  each  specimen  using  Micro-Measurement  WA-00-120WT-350 
strain  gages.  Young's  moduli  and  Poisson's  ratios  were  determined. 

Static  strengths  were  obtained  for  both  [Og]  and  [90g]  specimens 
in  the  wet  and  dry  conditions.  However,  a  series  of  strengths  were 
obtained  from  each  [90.J  specimen.  This  was  accomplished  by  testing 


TABLE  1.  SPECIMENS  TESTED 
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the  remaining  sections  after  each  failure.  Using  a  minimum  specimen 
length  requirement  of  1.5  inches,  up  to  eight  values  of  tensile 
strength  could  be  measured  for  each  initially  unbroken  12  inch  long 
[90g]  specimen.  Testing  the  [90g]  specimens  in  this  manner  allows  the 
experimentalist  to  isolate  the  first  or  "weakest  link"  strength  from 
the  subsequently  higher  strength  values.  Correlation  of  lamina  strength 
values  with  laminate  crack  spacing  equilibrium  models  developed  by 
Reifsnider  [17]  was  the  objective  of  the  weakest  link  [90g]  tests. 

Elastic  modulus  was  measured  ior  all  [90g]  tests  using  a  biaxial  strain 
gage  unless  the  previous  break  location  was  at  or  near  the  strain  gage. 

2.3  Preliminary  Tension  Tests  of  Qiiasi-Isotropic  Laminates 

The  objectives  of  these  static  tests  are  to  initially  isolate 

\ 

unique  free  edge  damage  states  in  type  I  and  type  II  laminates  and  show 
how  these  unique  damage  states  develop  into  the  characteristic  damage 
state  (CDS)  when  the  laminate  load  is  increased.  The  replica  technique 
as  utilized  by  Stalnaker  and  Stinchcomb  [13]  is  best  suited  for 
recording  and  observing  damage  which  develops  along  the  laminate  free 
edge.  By  using  the  replica  technique  an  image  of  damage  from  a  polished 
laminate  free  edge  can  be  imprinted  onto  a  thin  acetate  strip  while  the 
load  is  held  constant.  The  image  of  free  edge  damage  is  transferred  onto 
the  acetate  strip  under  pressure  by  simultaneously  wetting  the  polished 
laminate  surface  and  surface  of  the  acetate  with  acetone.  By  using  this 
technique  the  development  of  damage  can  be  recorded  by  taking  replicas  in 
increments  of  increasing  load.  One  edge  on  each  specimen  was  polished 
with  3  micron  alumina  oxide  particles  using  a  standard  metal lographic 


polishing  wheel  with  felt  cloth. 

The  specimens  used  in  the  preliminary  static  tests  listed  in  Table 
2  were  fabricated  by  Soutwest  Research  Institute.  Each  specimen  was 
incrementally  loaded  in  tension  by  a  hydraulic-load  controlled  Tinius 
Olsen  machine.  The  replica  technique  was  used  to  record  the  damage 
state  before  and  after  the  wet  and  dry  conditioning  as  well  as  at  the 
various  static  load  levels  listed  in  Table  2.  The  interpretation  of 
damage  imprinted  on  the  replicas  was  used  as  a  basis  for  improved  test 
procedures  used  in  the  quasi-isotropic  static  tension  tests  listed  in 
Table  3.  Strains  along  tlie  load  axis  were  measured  at  the  center  of 
each  specimen  using  Micro-Measurements  CEA-13-062VW-120  strain  gages. 

The  failure  strength  and  Young's  modulus  were  recorded  for  each  specimen 
which  was  incrementally  loaded  to  failure. 

2.4  Quasi-Isotropic  Static  Tension  Tests 

As  shown  in  Table  3  the  specimens  were  divided  into  two  groups. 

The  first  group  of  specimens  were  quasi-statically  loaded  to  failure 
with  no  replicas.  The  second  group  of  specimens  were  incrementally 
loaded  to  failure  with  replicas  taken  at  100  lb.  intervals.  The 
objective  of  the  quasi-statically  loaded  tests  was  to  statistically 
demonstrate  the  differences  in  wet  or  dry  laminate  strengths  for  type 
I  and  type  II  laminates.  As  discussed  in  the  introduction,  the  difference 
between  type  I  and  type  II  experimentally  determined  laminate  strengths 
provides  a  basis  for  investigating  how  the  free  edge  damage  develops 
under  load  and  influences  the  laminate  strength.  It  follows  that  the 
objective  of  the  second  group  of  tests  is  to  demonstrate  how  the 


TABLE  2.  PRELIMINARY  QUASI-ISOTROPIC  TENSION  TESTS 
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TABLE  3.  QUASI-ISOTROPIC  STATIC  TENSION  TEST 
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damage  in  wet  or  dry  quasi-isotropic  laminates  develop  into  the 
characteristic  damage  state  (CDS)  prior  to  fracture  of  the  laminate. 

In  order  to  characterize  the  development  of  damage,  replicas  at  100  lb. 
load  increments  were  obtained  for  only  one  specirsen  which  was  conditioned 
into  a  wet  or  dry  state  of  stress.  The  four  specimens  which  were  incre¬ 
mentally  loaded  to  failure  are  instrumented  with  biaxial  strain  gages 
located  at  the  center  of  each  specimen. 

2.5  Quasi-Isotropic  Fatigue  Tests 

The  objective  of  the  fatigue  tests  is  to  characterize  the  develop¬ 
ment  of  damage  in  wet  and  dry  quasi-isotropic  laminates  when  cyclically 
loaded.  As  shown  in  Table  4,  replicas  were  taken  for  two  types  of  cyclic 
loads.  The  coaxing  cyclic  load  when  used  with  replicas  implies  that 
the  experimentalist  can  interact  with  the  fatigue  test  such  that  a 
unique  state  of  damage  could  be  isolated  or  "coaxed"  out  of  the  laminate. 
When  the  steady  state  cyclic  load  is  used  the  experimentalist  records 
the  regular  sequence  of  damage  events  on  replicas  and  does  not  interact 
with  the  fatigue  test.  The  remaining  fatigue  tests  are  not  replicated, 
but  form  a  sufficient  data  base  for  studying  the  overall  laminate 
response  (i.e.  residual  strength,  stiffness  change)  when  using  either 
the  coaxing  or  steady  state  cyclic  loads.  Initial  Young's  modulus  is 
measured  using  a  MTS  clip-on  extensometer  while  the  specimen  is  being 
quasi-statically  preleaded  up  to  the  maximum  cyclic  load.  All  tests 
were  cycled  at  10  Hertz  on  a  load  controlled  MTS  machine.  For  both 
coaxing  and  steady  state  cyclic  loads  the  change  in  Young's  moduli 
is  periodically  measured  by  stopping  the  test  and  reloading  up  to 
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the  maximum  cyclic  load.  A  constant  stress  ratio  of  0.1  v/as  used  for 
all  fatigue  tests.  The  drying  out  of  all  wet  specimens  during  the 
fatigue  tests  was  eliminated  by  simply  wrapping  a  wet  paper  towel 
around  the  test  specimen. 

In  order  to  coax  the  desired  unique  free-edge  damage  state  the 
experimentalist  should  be  able  to  stop  the  test  momentarily  and  observe 
the  development  of  damage.  The  author  discovered  that  the  image  of 
damage  replicated  onto  an  acetate  strip  can  be  immediately  viewed  on 
a  microfiche  card  reader  with  sufficient  magnification  and  clarity 
for  interpretation.  The  experimentalist  can  then  interpret  the  state 
of  damage  during  the  test  and  decide  which  test  variables  should  be 
changed  in  order  to  furthe*"  isolate  or  alter  the  observed  damage.  As 
a  result  of  this  interactive  test  technique,  no  clearly  defined  test 
procedure  is  established.  Although  the  test  procedure  is  open  ended, 
the  experimentalist  must  base  his  interpretation  of  damage  on  a  rational 
mechanistic  approach.  This  requires  an  understanding  of  the  initial 
damage  free  stress  state  and  the  CDS  which  will  be  explained  in  Chapter 
III.  Therefore  the  guidelines  or  rationale  used  to  interpret  the  damage 
state  when  using  the  coaxing  cyclic  load  is  outlined  in  Section  4.4.1 
following  Chapter  III.  One  specimen  from  each  wet  and  dry  stress  state 
was  fatigued  by  a  coaxing  cyclic  load.  Replicas  at  200  lb.  intervals 
and  residual  strengths  were  obtained  by  incrementally  loading  the 
specimens  to  failure. 

The  objective  of  the  steady  state  cyclic  load  tests  is  to  record 
the  regular  sequence  of  damage  events  which  develop  into  the 
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characteristic  damage  state.  Replicas  are  taken  at  200  lb.  intervals 
during  the  initial  preload  up  to  the  maximum  loads  listed  in  Table  4. 
Replicas  were  also  taken  at  5K,  lOK,  50K,  and  lOOK  cycles.  For  each 
wet  and  dry  stress  state,  two  specimens  were  cycled  to  lOOK  cycles  and 
one  specimen  cycled  to  IM  cycles.  Residual  strengths  were  obtained 
after  lOOK  or  IM  cycles  with  replicas  taken  at  100  lb.  intervals. 


III.  STRESS  ANALYSIS 


3.1  Linear  Elastic  Homogeneous  Transversely  Isotropic  (LEHTI) 

Material  Properties 

For  fiber-reinforced  composite  laminates  it  has  been  experimentally 
demonstrated  [13]  that  damage  initiates  near  the  laminate  free  edge 
in  the  form  of  ply  cracking  and  interply  delaminations.  Finite  element 
and  finite  difference  models  and  perturbation  solutions  [34,24,25]  wer.* 
developed  in  order  to  predict  stresses  near  the  free  edge  prior  to  the 
formation  of  damage.  Both  the  finite  element  and  finite  difference 
models  require  all  nine  orthotropic  lamina  properties.  If  the  initial 
damage-free  stress  state  existing  near  the  free  edge  uniquely  influences 
the  subsequent  formation  of  damage  then  accurate  out-of-plane  elastic 
properties  (e.g.  v2-j  and  G22)  must  be  used  when  evaluating  this  stress 
state.  Unfortunately  all  nine  elastic  lamina  prope.* ties,  are  difficult 
to  obtain;  therefore  it  is  not  unconEnon  to  make  simplifying  approxima¬ 
tions,  such  as  ^23”'’’12”^13  ^23~^12'^13’  order  to  obtain  solutions. 

A  parametric  evaluation  by  Kriz  [26]  of  these  approximations  and  other 
property  variations  on  the  interlaminar  stresses  in  angle-ply  laminates 
demonstrates  that  more  emphasis  should  be  placed  on  accurately  calcu¬ 
lating  or  experimentally  measuring  all  lamina  elastic  properties  if 
meaningful  stress  distribtuions  are  to  be  obtained  from  the  various  models. 

In  this  investigation  each  laminate  layer  (lamina)  is  assumed  to  be 
constructed  of  fiber  reinforced  materials  as  shown  in  Figure  1.  When 
the  fibrous  lamina  are  assumed  to  be  homogeneous  and  transversely 
isotropic,  the  number  of  independent  material  constants  reduces  to 


20 


21 


FIGURE  1.  COORDINATE  SYSTEM  DEFINITION 
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five.  The  most  commonly  used  engineering  constants  are  listed  below: 

E,.  E2  =  E3,  G,2  =  G^3,  =  v^3,  633 

The  symbols,  E,  G,  and  v  are  Young's  modulus,  shear  modulus,  and 
Poisson's  ratio,  respectively,  and  numbered  subscripts  refer  to  the 
lamina  coordinates,  shown  in  Figure  1.  For  transversely  isotropic 
materials,  the  Poisson's  ratio  in  the  2-3  plane  is  related  to  E2 
and  G23  by  Equation  (1) 

V23  =  E2/2G23-I  (1) 

Reasonable  estimates  for  ^23  and  G23  together  with  other  laminate 
elastic  properties  have  been  obtained  experimentally  using  mechanical 
and  ultrasonic  techniques  [18,19,20].  A  complete  set  of  elastic 
lamina  properties,  including  v^3  and  G23.  can  be  calculated  using 
equations  developed  by  Hashin  [27]  as  listed  in  Appendix  B  given  the 
isotropic  matrix  and  transversely  isotropic  fiber  elastic  properties. 
Unfortunately  the  complete  set  of  transversely  isotropic  graphite 
fiber  properties  are  difficult  to  obtain  due  to  the  small  fiber 
diameter. 

Oean  and  Turner  [18]  demonstrated  that  most  of  the  graphite  fiber 
properties  could  be  determined  by  curve  fitting  ultrasonically 
evaluated  stiffnesses  over  a  range  of  fiber  volume  fraction  by  using 
Hashin's  equations  v/ritten  in  terms  of  ultrasonic  stiffnesses. 
Unfortunately,  due  to  scatter  in  ultrasonic  data.  Dean  and  Turner 
were  unable  to  extrapolate  for  all  fiber  properties.  Kriz  and 
Stinchcomb  [20]  improved  on  the  extrapolation  technique  of  Dean  and 
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Turner  and  showed  that  although  ultrasonic  data  are  scattered,  all 
fiber  properties  can  be  determined  and  that  five  independently  extra¬ 
polated  graphite  fiber  properties  are  consistent  with  the  assumption 
that  fibers  can  be  modeled  as  transversely  isotropic.  However,  there 
is  an  error  in  the  equations  used  by  Kriz  and  Stinchcomb  [20]  to 
calculate  values  of  622.  The  correct  expression  for  62  is  shown 
below  in  Equation  (2)  where  62  is  used  in  Equation  B.2(b)  of 
Appendix  B  to  calculate  G22"  Fortunately  all  graphite/epoxy 

Kf  ^ 

®2  +^2G7~ 

^tt  -tt 

lamina  properties  whicr.  depend  on  the  correction  for  62  are  changed  by 
less  than  1  percent.  It  is  also  noted  that  the  corrections  have  no 
effect  on  the  extrapolated  fiber  properties.  ■  Therefore,  all  observa¬ 
tions  and  conclusions  by  Kriz  and  Stinchcomb  [20]  are  unaltered.  Although 
these  corrections  are  not  numerically  significant  for  the  graphite/epoxy 
materials  evaluated  in  Ref.  [20],  there  could  be  other  fiber/matrix 
systems  for  which  these  corrections  could  becom.e  significant. 

Graphite/epoxy  lamina  properties  were  also  experimentally  measured 
by  Ishikawa  et  al  [19]  using  mechanical  test  techniques  and  fiber 
properties  v^ere  extrapolated.  Comparison  of  extrapolated  graphite 
fiber  properties  along  with  epoxy  matrix  properties  is  demonstrated 
in  Table  5.  A  final  sot  of  fiber  and  matrix  elastic  properties  are 
chosen  and  listed  in  Table  6. 

Since  fibers  properties  are  assumed  unchanged  by  moisture. 
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TABLE  6.  WET  AND  DRY  ELASTIC  PROPERTIES 
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the  wet/dry  lamina  properties  listed  in  Table  6  are  calculated  by 
substituting  the  fiber  properties  together  with  the  wet/dry  matrix 
properties  into  the  equations  listed  in  Appendix  B.  The  changes  in 
elastic  matrix  properties  due  to  moisture  absorbtion  are  chosen  such 
that  the  corresponding  changes  in  elastic  lamina  properties,  reflect 
experimental  differences  in  wet/dry  lamina  properties.  For  example, 
an  11  percent  reduction  in  lamina  modulus,  E2,  is  calculated  when  a 
25  percent  decrease  in  matrix  modulus,  E^,  together  with  a  27  percent 
reduction  in  matrix  shear  modulus,  G^,  are  substituted  into  the  equa¬ 
tions  of  Appendix  B.  The  predicted  decrease  of  11  percent  for  E2  is 
comparable  with  experimental  reductions  of  10.7  percent  by  Hahn  [28] 
and  11.2  percent  by  Hofer  et  al  [29].  A  decrease  of  0.5  percent  is 
predicted  for  E-j  v/hich  does  not  compare  well  with  a  decrease  of  4.5 
percent  measured  by  Hahn  and  Kim  [28]  and  16.3  percent  measured  by  Hofer 
[29].  An  insufficient  data  base  on  wet/dry  T300/5208  lamina  shear 
properties  does  not  provide  meaningful  comparison  with  the  predicted 
reductions  of  23  percent  for  G^2  percent  for  62^.  To  date  no 

reliable  and  reproducible  shear  test  method  has  been  widely  accepted 
to  verify  the  trends  predicted  for  G12  and  G22. 

Variations  in  test  techniques  together  with  variations  in  test 
specimens  due  to  manufacturing  processes,  quality  control,  etc.  are 
obvious  problems  which  account  for  much  of  the  variation  in  experimental 
data.  Variations  in  fiber  volume  fraction,  temperature  and  duration  of 
cure  [22]  could  provide  additional  explanations  for  variations  in 
experimental  data.  For  these  reasons  the  author  has  choosen  to  select 


lamina  properties  as  calculated  from  the  equations  in  Appendix  B  which 
model  the  change  in  experimental  data  due  to  nxjisture  absorption  as  a 
reduction  of  matrix  elastic  properties.  The  experimental  value  of 
18.8  Msi  was  chosen  for  corresponding  to  a  fiber  volume  fraction  of 
55  percent  as  calculated  from  equations  in  Appendix  B  when  fiber  and 
wet  matrix  properties  were  chosen.  Wet  and  dry  type  I  and  type  II 
laminate  stiffnesses  were  calculated  from  laminate  olate  theory  using 
the  wet  and  dry  lamina  properties  at  55  percent  fiber  volume  fraction. 
To  be  consistent  with  pi-evious  analysis  [16]  a  value  oT  0.65  msi  was 
chosen  for  the  shear  transfer  layer  modulus,  Interestingly,  a 

value  of  0.65  msi  is  calculated  for  G-i^  v/hen  at  a  fiber  volume  fraction 
of  42  percent  is  used  in  Equation  (13)  of  Appendix  B.  It  is  reasonable 
to  assume  that  the  shear  transfer  which  acts  over  a  small  distance 
above  and  below  a  resin  rich  interface  should  have  a  slightly  lower 
fiber  volume  fraction  when  compared  to  55  percent  in  the  plies  adjacent 
to  the  interface. 

The  fiber,  matrix,  lamina,  laminate,  and  interface  elastic 
properties  were  chosen  to  demonstrate  how  variations  in  wet  and  dry 
conditioning  of  type  I  and  type  II  laminates  affect  the  free  edge 
stresses  and  subsequent  development  of  dam=*ge  at  the  free  edge. 
Variations  in  experimentally  measured  elastic  lamina  properties  can 
be  large.  For  this  reason  the  lamina  elastic  properties  were  chosen 
such  that  only  variations  due  to  moisture  absorption  are  modeled 
using  the  equations  of  Appendix  B.  Elastic  properties  chosen  in 
this  way  do  not  include  the  inherent  variations  in  experimentally 


neasured  elastic  properties.  These  properties  provide  a  consistent  set 
of  elastic  properties  to  be  used  in  the  2-D,  3-D,  and  crack  spacing 
models  following  this  section. 

3.2  Two  Dimensional  Thin  Laminate  Theory  Stress  State 

The  first  approximation  of  a  laminate  stress  state,  neglecting 
out-of-plane  stresses  (plane  stress  is  assumed),  is  obtained  by 
laminated  plate  theory.  Since  the  laminate  thickness  is  much  smaller 
than  its  other  dimensions  Kirchhoff's  hypothesis  can  be  used.  As  a 
result  of  these  approximations  each  layer  of  the  laminate  is  assumed 
to  be  in  a  state  of  plane  stress  v/hich  can  be  directly  related  to 
inplane  loads  and  moments.  A  complete  development  of  the  classical 
thin  laminated  plate  theory  is  given  by  Jones  [30].  The  objective  of 
this  section  is  to  predict  the  stresses  existing  in  the  thin  laminate 
plane  prior  to  the  damage  event,  using  laminate  theory,  aru'  assuming 
no  initial  damage  exists.  The  predicted  stresses  are  the  result  of 
an  applied  in-plane  mechanical  load  acting  together  with  residual 
stresses  due  to  swelling  from  moisture  absorption  and  curing  at 
elevated  temperatures. 

The  combined  effect  of  a  mechanical  load,  N^,  and  residual  stresses 
R 

0^,  on  the  individual  ply  stresses,  o.,  is  shown  in  Equation  (3)  using 
the  notation  of  Kim  and  Hahn  [11]. 

Of  =  CijAjjN,  (3) 

where 

Q^j  =  reduced  ply  stiffnesses 
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* 


A..  =  laminate  stiffnesses 
^  *3 

N|^  =  laminate  stress  resultants 

p 

0^  =  residual  ply  stresses 

p 

The  residual  ply  stresses,  0!/,  can  be  represented  in  terms  of 
differfences  in  lamina  and  laminate  thermal  and  moisture  strains. 


oM  M, 


(4) 


where 

=  laminate  thermal  strain 
3 

=  laminate  moisture  strain 
3 

el  =  lamina  thermal  strain 
3 

M  •  • 

e.  =  lamina  moisture  strain 

These  strains  can  be  defined  in  terms  of  coefficients  of  expansion. 


where 

а.  =  lamina  thermal  coefficient  of  expansion 
3 

б.  =  lamina  moisture  coefficient  of  expansion 
3 

AT  =  change  in  temperature 

aM  =  percent  weight  gained  by  moisture  absorption 

a.  =  laminate  thermal  coefficient  of  expansion 
3 

6-  =  laminate  moisture  coefficient  of  expansion 
3 

Laminate  coefficients  of  expansion  are  defined  by  integrating  the 
lamina  strains  written  in  terms  of  coefficients  of  expansion  through 
the  laminate  thickness,  {-h/2,  h/2),  as  shown  below. 
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/h/2 

-h/2 


(6) 


/  h/2 
-h/2 


For  the  type  I  and  type  !I  laminates,  damage  first  occurs  in  the 
90"'  plies  in  the  form  of  cracks  transverse  to  the  direction  of  the 
uniaxial  laminate  tension  load,  N^.  Since  cracks  first  occur  in  the 
90*^  plies  we  calculate  the  stresses  in  this  ply  by  expanding  equation 
(3) ,  as  shown  below. 


Ov  P  h  +  a, 

*  A  A  ^  ' 

^xx  yy  ^xy 


•Q-i/sA  p. 

^  ^  ^22  yy  ^12  xy  „  .  R 

Q  _  - ^ ^  ^  ^ 


A  A  -A" 
XX  yy  xy 


(8) 

(9) 


0g  =  0  (10) 

where  subscripts  1  and  2  refer  to  the  lamina  coordinates  defined  in 
Figure  1  and  subscript  6  is  the  contracted  notation  for  the  lamina 
shear  stress  in  the  1-2  plane.  The  initial  crack  in  the  90°  ply  is 
labeled  FPF  for  first  ply  failure  and  the  corresponding  laminate  load 
is  labeled,  .  Although  0^  and  0^  both  exist  in  the  90°  ply  at  FPF 
the  transverse  stress,  02,  has  reached  the  transverse  strength  while 
is  much  lower  than  the  longitudinal  strength.  Therefore  a  maximum 
stress  failure  criterion  can  be  used  which  assumes  FPF  occurs  when  the 
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value  of  02  in  equation  (9)  is  equal  to  the  unconstrained  uniaxial 

tensile  strength,  T,  of  a  90°  test  specimen.  Using  the  maximum  stress 

FPF 

failure  criterion,  equation  (9),  can  be  used  to  solve  for  by 


replacing  02  with  T. 


FPF  ^  (T-c^) 

X  Q22\y"^12\y  ^ 


(11) 


Prediction  of  the  FPF  laminate  load,  depends  on  how  accurately 

R 

02  is  predicted.  From  equation  (5)  the  value  calculated  for  02  depends 
upon  approximations  made  for  lamina  thermal  and  moisture  coefficients 
together  with  realistic  estimates  of  percent  weight  gained  due  to 
moisture  absorption,  AM,  and  changes  in  temperature,  aT,  from  the  stress 
free  temperatute.  By  measuring  the  warpage  of  unsymmetric  [±e] 
laminates  Pagano  and  Hahn  [31]  suggest  a  stress  free  cure  temperature 
of  250°F  for  T300/5208  which  is  considerably  lower  than  the  350°F  cure 
temperature.  For  T300/5208  graphite/epoxy  Hahn  [^8]  measured  lamina 
thermal  and  moisture  coefficients  of  expansion,  as  shown  below. 


a.|  =  -0.17  ne/°F;  02  =  15.6  yE/°F,  a  0;  ^2  =  ^900  v»c/^!;H20 

(12) 

where  ^2  is  calculated  from  experimental  data  assuming  a  moisture 
threshold,  C^,  of  0.4%  as  shown  in  Figure  2.  Independently  Crossman 
[10]  demonstrated  that  swelling  of  T300/5208  graphite/epoxy  exhibits 
the  same  moisture  threshold  with  less  scatter  of  data  and  measured 
lamina  thermal  and  moisture  coefficients  of  expansion  shown  belcw. 


=  0.16  pe/°F;  02  "  14.3  uc/°F;  p.j  ~  0‘,  $2  ~  5000  u£/'=H20  (13) 
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The  coefficients  of  thermal  expansion  shown  below  were  the  same  as 
used  by  Nagarkar  and  Herakovich  [32]  and  moisture  expansion  coefficients 
were  obtained  by  curve  fitting  experimental  data  reported  by  Crossman 
[10]  assuming  a  moisture  threshold  of  0.4%. 

=  -0.23  ue/°F;  =  14.9  viE/°F;  e  0;  6^  =  5560  pe/XH^O 

(14) 

The  laminate  thermal  and  moisture  coefficients  of  expansion  are 
calculated  for  type  I  and  type  II  laminates  by  substituting  the  wet  and 
dry  elastic  properties  listed  in  Table  6  together  with  the  lamina 
expansion  coefficients  (14)  into  equations  (6)  and  (7),. 

“XDRY  "  “YDRY  “  ve/°F;  =  484  yc/%H20 

. “XWET  ""  “YWET  ^^XWET  ^'WET  '  ut7..H20 

It  is  now  possible  to  calculate  wet  and  dry  residual  stresses  at 
room  temperature  by  using  equation  (5)  where  AM  was  measured  as  1.2 
pe)'cent  in  the  fully  saturated  state  and  a  temperature  change  of  -180°F 
was  estimated  using  Hahn's  approximation  of  250'F  for  the  stress  free 
temperature.  The  Q.j  were  calculated  using  wet  and  dry  properties 
1 isted  in  Table  6. 

*^2DRY  ~  3.43  ksi  ^2WET  ~  ~1‘^5  ksi  (15) 

The  laminate  FPF  load  is  now  calculated  from  equation  (9)  using 
elastic  properties  in  Table  6  together  with  the  difference  between  the 
transverse  strength,  T,  and  the  residual  stress  state.  The  values  for- 
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the  weakest  or  first  transverse  strengths,  T,  as  described  in  Section 

FPF 

2.2  are  listed  in  Table  10  and  used  in  Equation  (11)  to  calculate  . 

Using  only  dry  properties  in  equation  (11),  a  dry  FPF  laminate  load, 
^XDKY’  calculated.  Similarly  a  wet  FPF  laminate  load,  was 

calculated  using  only  wet  properties. 

'^XDRY  "  ^^WET  " 

The  FPF  laminate  loads  for  type  I  and  type  II  laminates,  predicted  by 
laminate  theory,  are  identical  since  changing  the  stacking  sequence 
has  no  effect  on  o^. 

It  is  interesting  to  note  that  the  predicted  wet  residual  stress 
state  02y£y  in  the  90°  ply  is  in  compression  by  -1.96  ksi.  For  this 
reason  a  higher  laminate  load  is  required  to  produce  FPF  when  the 
laminate  has  absorbed  1.2  percent  moisture  by  weight  gain. 

Although  this  section  vjas  primarily  concerned  with  evaluating  the 
stress  state  in  the  90°  ply  prior  to  FPF,  all  0°  and  45°  ply  stresses 
in  the  wet  or  dry  state  have  also  been  calculated  using  equation  (3) 

ppp  ppp 

in  similar  fashion.  The  calculated  ply  stresses  for  ^^XWET 

are  listed  in  Tables  7  and  8  for  comparison  with  the  predicted  finite 
element  stresses  in  Section  3.3  which  use  the  sane  wet  and  dry  FPF 
laminate  loads. 

3.3  Three  Dimensional  Finite  Element  Theory  (LEHTI)  Stress  State 

The  objective  of  this  section  is  to  predict,  using  finite  element 
theory,  the  stresses  existing  in  the  interior  and  along  the  edge  of 
type  I  and  type  II  laminates  prior  to  the  first  carnage  event.  The 


predicted  stresses  are  the  result  of  a  wet  or  dry  FPF  laminate  load, 
which  were  calculated  in  section  3.2,  acting  together  with  residual 
stresses  due  to  swelling  from  moisture  absorption  and  shrinking  after 
cooling  down  from  an  elevated  stress  free  temperature. 

The  Finite  Element  Model  (FEM)  used  in  this  section  was  developed 
at  Virginia  Tech  in  the  Engineering  Science  and  Mechanics  Department 
by  graduate  students  under  the  direction  of  Dr.  Carl  Herakovich.  Several 
thesis  and  dissertations  [32,33,34,35]  have  been  involved  in  the  develop¬ 
ment  of  this  FEM.  The  current  version  of  the  FEM  {NONCOM  III)  was  used 
in  this  investigation.  The  general  formulation  of  this  FEM  is  reported 
in  the  most  recent  reference  [32]. 

The  FEM  represents  the  laminate  in  a  state  of  generalized  plane 
strain  in  the  x  load  direction  by  using  constant  strain  triangles  as 
shown  in  Figure  3.  Because  of  symmetry  conditions  for  type  I  and  type 
II  laminates  only  a  quarter  of  the  y-z  plane  is  modeled.  Boundary 
conditions  are  imposed  such  that  all  nodes  along  y  =  0  are 
constrained  from  moving  in  the  y-direction  but  are  free  to  move  in 
the  z-direction,  and  all  nodes  along  z  -  0  are  constrained  from 
moving  in  the  z-direction  but  are  free  to  move  in  the  y-direction. 

The  coninon  node  at  y  =  0,  z  =  0  is  held  fixed.  The  externally 
applied  nodal  forces  along  Z  =  H  and  y  =  B  are  prescribed  to  be  zero. 
These  prescribed  nodal  forces  represent  fie  free  edge  and  free 
surface  stress-free  boundary  conditions.  The  FEM  grid  used  to  model 
the  four  layer  quarter  plane  as  shown  in  Figure  3  uses  763  elements 
and  438  nodes  with  96  triangular  elements  at  the  laminate  edge.  As 
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shown  in  Figure  3  the  size  of  the  elements  are  reduced  near  the  free 
edge  in  order  to  model  the  stress  gradients  in  this  region.  The 
smallest  edge  element  used  in  this  dissertation  is  shown  in  Figure 
3  such  that  15-20  fibers  fill  the  triangular  area  with  at  least  3 
fibers  along  the  triangular  edge. 

Although  the  finite  element  method  is  not  limited  by  how  small 
the  elements  are  chosen,  the  smallest  element  was  chosen,  in  this 
investigation,  such  that  the  composite  material  modeled  within  the 
smallest  element  shown  in  Figure  3  can  be  assumed  to  behave  as  a 
homogeneous  material.  This  is  only  a  limitation  on  the  idealization 
of  homogeneity  used  in  the  constitutive  relations  by  the  finite 
element  method. 

The  author  chose  not  to  model  the  nonlinearities  of  T300/5208 
since  Crossman  [10]  experimentally  verified  that  the  calculated 
linear  elastic  residual  stress  state  can  be  used  as  an  approximation 
for  predicting  the  curvature  of  unsymmetric  T300/5208  laminates. 
Although  the  differences  between  the  experimental  and  predicted 
curvatures  are  most  likely  due  to  nonlinearities  of  T300/5208,  these 
differences  are  small  and  are  ignored  in  order  to  minimize  the  number 
of  approximations  introduced  into  an  idealized  FEFi.  The  point  is 
made  that  these  nonlinearities  are  not  insignificant-;,  but  that  there 
is  much  yet  to  be  explained  about  FPF,  as  will  be  demonstrated'  in 
Section  4.5,  when  using  only  the  linear  elastic  idealization.  A 
case  in  point  is  the  crack  spacing  model  developed  by  Reifsnider 
et  al  [16]  which  uses  an  idealized  linear  elastic  shear  lag  model 
to  predict  the  crack  spacing  following  the  FPF  event. 


Only  the  linear  elastic  facility  of  the  FEM  (NONCOM  III)  is 
utilized  by  using  the  wet  and  dry  elastic  properties  listed  in  Table  6 
together  with  the  lamina  expansion  coefficients  (14)  from  section  3.2 
where  AT  -  -180°F  and  AM  =  1.2  percent.  Since  the  present  FEM  is  not 
formulated  with  respect  to  Hahn's  moisture  threshold,  c^,  the  ^2 
expansion  coefficient  is  recalculated  to  give  the  same  maximum  swelling 
strain  as  shown  in  Figure  2.  For  comparison  with  the  FPF  stress  state 
calculated  in  section  3.2,  the  same  wet  and  dry  FPF  laminate  loads 
are  used  by  the  FEM  to  predict  the  edge  stress  state.  In  order  to 
study  the  combined  effect  of  the  moisture,  temperature,  and  FPF 
laminate  load  on  the  edge  stresses,  the  FEM  calculates  the  edge 
stresses  due  to  moisture,  temperature,  and  FPF  load  separately  and 
then  superposes  the  stresses  to  idealize  the  final  wet  and  dry 
edge  stress  state  for  type  I  and  type  II  laminates. 

Although  it  may  be  instructive  to  demonstrate  how  the  initial 
residual  dry  stresses  changes  with  stacking  sequence,  followed  by 
moisture  absorption  and  mechanical  load,  only  the  final  superposed 
wet  or  dry  edge  and  interior  stress  states  existing  prior  to  damage 
formation  are  presented  in  the  section.  Crossman  [10,14]  has  already 
provided  an  interesting  evaluation  on  how  various  stacking  sequences, 
uniformily  distributed  moisture  levels  and  applied  loads  affect  the 
edge  stress  state.  The  emphasis  in  the  present  investigation  is  to 
compare  the  final  superposed  wet  or  dry  edge  and  interior  stress 
states  existing  prior  to  damage  formation  with  the  first  formation 
of  damage  as  observed  using  the  replica  technique. 
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The  free  edge  and  interior  stress  state  for  a  type  I  wet  laminate 

with  a  wet  FPF  laminate  load  is  graphically  demonstrated  in  Figures 

4-11.  A  four  element  averaging  scheme  outlined  in  reference  [35]  is 

used  to  calculate  average  and  stress  distributions  along 

the  laminate  interfaces  as  shewn  in  Figures  4-6.  The  same  four 

element  averaging  scheme  is  used  to  average  o  ,  o  ,  and  t  stresses 

X  y  xy 

along  the  midlayers  as  shov/n  in  Figures  7-9.  Through  the  thickness 
distributions  plotted  in  Figures  10  and  11  use  a  two  element  averaoing 
scheme. 

Only  stress  distributions  for  type  I  wet  are  shown  in  Figures 

4-il.  The  edge  and  interior  values  of  stress  for  the  type  I  dry, 

type  II  dry  and  type  II  wet  laminates  af^  summarized  in  Tables  7  and 

8  along  with  the  laminate  plate  theory  stress  state  for  comparison. 

As  a  check  all  interior  FEM  stresses  are  eqiivalent  to  stresses 

predicted  by  laminate  theory,  as  shown  in  Tables  7  and  8.  Minor 

variations  in  the  interior  stresses  between  FEM  and  laminate  theory 

exist  since  laminate  theory  models  the  o,  ,  a  ,  and  t  stresses  as 

X  y  xy 

constants  v/ithin  the  layer  v/ith  a  discontinuity  in  stress  between 
layers. 

For  later  reference  with  respect  to  formation  of  edge  damage, 
the  damage  free  and  stress  state  for  v/et  and  dry  type  I  and 
type  II  lamiiates  are  compared  in  Figure  12. 

3.4  Characteristic  Damage  State  Model  Definition 

The  characteristic  damage  state  (CDS)  was  first  introducted  by 
Reifsnider  et  al  [16]  as  a  laminate  property  which  could  be  defined 


43 


44 


* 


Y/B 


FIGURE  6.  VARIATION  OF  NEAR  THE  FREE  EDGE  OF  A  TYPE  I  WET 
LAMINATE  (N^=1900  Ib/in,  .■.T=-180"F,  ) 
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FIGURE  7.  VARIATION  OF  NEAR  THE  FREE  EDGE  OF  A  TYPE  I  WET 
LAMINATE  (N^=I900  Ib/in,  &T=-1S0°F,  aM=1.250 


4. 


o  c  o  Stresses  (ksi) 
X  y  z 


FIGURE  10.  VARIATION  OF  7^  AND  THROUGH  THE  THICKNESS  AT 
Y/B=0.998  FOR  A  TYPE  I  WET  LAMINATE  (N^=190C  Ih/in. 
aT=-18C'’F,  1.2'‘>) 


r 
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FIGURE  11.  VARIATION  OF  THROUGH  THE  THICKNESS  AT 

Y/B=0.998  FOR  A  TYPE  I  WET  LAMir;ATE  (N^=1900  Ib/in, 
at=-180‘^f,  An=i.2':.) 
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FIGURE  12.  COMPARISON  OF  AND  FPF  STRESSES  IN  WET  AND  DRY 
TYPE  I  AND  TYPE  II  LAMINATES  AT  Y/3=0.998 


as  a  predictable  and  regular  pattern  of  ply  cracks  which  exists  prior 
to  fracture.  Definition  of  the  damage  state  prior  to  fracture  is 
facilitated  by  a  shear  lag  model,  employed  in  reference  [16],  in 
order  to  predict  the  stable  patterns  of  ply  cracks  which  develop  after 
the  FPF  event.  Although  other  damage  in  the  form  of  interply  cracks 
and  longitudinal  ply  cracks  exists,  the  stable  pattern  of  ply  cracks 
form  the  basis  of  the  CDS.  Understanding  how  this  crack  pattern  and 
subsequent  stresses  are  distributed  in  this  CDS  prior  to  fracture 
should  explain  the  difference  in  fracture  strengths  for  type  I  and 
type  II  laminates.  Talug  [36]  investigated  the  states  of  stress 
surrounding  a  partially  damaged  r'^'^’cn  in  type  I  and  type  II  laminates 
and  demonstrated  the  differences  that  exist  in  the  response  of  these 
laminates  even  in  regions  remote  from  the  edges. 

The  shear  lag  model  duveioped  in  references  [37,33]  predict  crack 
spacing  in  90°  plies  with  increasing  load.  This  shear  lag  model  was 
developed  for  the  specific  case  where  the  plies  adjacent  to  the  90° 
plies  are  of  the  same  orientation.  Since  this  restriction  prohibits 
the  use  of  these  models  for  a  type  II  laminate,  where  the  90°  plies 
are  constrianed  by  0°  and  45°  plies,  the  more  general  shear  lag  mode! 
formulated  by  Re’fsnider  et  a1  [16]  is  used  in  this  investigation. 

The  concept  of  a  characteristic  crack  pattern  which  is  introduced  by 
Reifsnider  et  al  [16]  is  also  utilized  in  this  investigation. 

Although  much  work  is  being  done  to  further  characterize  how 
the  CDS  uniquely  influences  the  laminate  response,  the  effort  in  this 
investigation  is  to  better  define  the  CDS  with  respect  to  the  effect  of 


moisture.  Since  the  stable  pattern  of  cracks  is  the  basis  of  the  CDS, 
the  effects  of  moisture  on  the  crack  patterns  will  be  demonstrated. 

The  shear  lag  model  is  developed  by  Reifsnider  et  al  [16]  in  sufficient 
detail  such  that  only  basic  principles  are  restated  below. 

The  shear  lag  model  assumes  that  when  a  FPF  crack  forms  in  an 
off  axis  ply  the  stress  in  the  ply  is  recovered  from  zero  at 
the  crack  surface  to  the  unbroken  ply  stress  as  shown  in  Figure  13. 

The  stress  near  the  FPF  is  no  longer  carried  through  the  cracked  ply 
but  is  transfered  by  shear  through  plies  adjacent  to  the  cracked  ply 
until  the  stress  in  the  cracked  ply  has  recovered  to  99.9  percent  of 
the  unbroken  stress  which  results  in  the  formation  of  a  second  crack. 
The  distance  between  the  first  and  second  crack  is  referred  to  as  the 
crack  spacing.  Cracks  continue  to  form  in  this  manner  with  increasing 
laminate  load  until  no  new  cracks  form  at  which  point  an  equilibrium 
crack  spacing  is  reached.  Ideally  the  cracked  ply  carries  only  99.9% 
of  the  unbroken  stress  until  equilibrium  crack  spacing  is  realized 
and  then  the  cracked  ply  no  longer  carries  any  stress.  The  crack 
spacing  is  said  to  be  characteristic  of  the  laminate  since  ply  and 
laniinate  elastic  properties  are  used  together  vn'th  the  ply  strength 
to  calculate  distance  between  cracks.  The  effect  of  the  stacking 
sequence  on  the  crack  spacing  is  included  in  the  shear  lag  model  by 
the  constraining  effect  of  adjacent  plies  on  the  redistribution  of 
stress  near  the  crack. 

For  type  I  laminates  the  redistributior,  of  stress  in  the  crack 
ply  due  to  a  shear  transfer  is  formulated  by  Reifsnider  et  al  [16] 
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in  more  general  terms  by  using  normalized  displacement,  u,  where  the 
governing  differential  equations  shown  below  are  the  result  of  an 
equilibrium  element  analysis  analogous  to  the  approach  used  in 
references  [39,40]. 


Boundary  conditions 


d  Uqp 

^  .2  *^45  ‘'so  ^ 

dx 

(17) 

d^U45 

®  .1^  ''so  V"  "  ^“45  ® 

dx 

(18) 

(X  .  »)  .  1 

(19) 

(20) 

(X  •  0)  .  0 

(21) 

U45  (x  -  0)  -  0 

(22) 

A  =  (b  Egg)/(2a  E^) 

(23) 

B  =  (b  c  E^5)/(a^  E^^; 

(24) 

U  =  [a^a/ /EyGj-j^lu 

(25) 

X  =  ^  ^ 

(26) 

where 


X 

=  normalized  distance  from  crack  surface 

A 

^90 

=  normalized  displacement  in  90°  ply 

“45 

=  normalized  displacement  in  45°  ply 

u 

=  displacement 

X 

=  distance 

=  stress  applied  to  laminate 

=  laminate  stiffness 

m 

<0 

0 

=  90°  ply  stiffness 

Sl 

=  shear  transfer  layer  shear  modulus 

a 

=  90°  ply  thickness 

b 

=  shear  transfer  layer  thickness 

c 

=  45°  ply  thickness 

Solutions  to  the  differential  equations  are  given  below, 

“90  ■  *^1^  ^ 

(27) 

U45  =  X  .  D^cie-  ^  D2C2e-^^ 

(28) 

where 

1/2 

rB+2A  1  fR^  + 

a,g  =  ±  ^  +  4A  1  J 

(29) 

2 

=  1  -  Aa 

(30) 

Cl  -  -  D2-/^s0i  -  002) 

(31) 

D„  -  1  -  Ab^ 

(3?) 

C2  ”  Di/lsO-j  -  C1D2) 


(33) 


A  similar  set  of  solutions  for  the  displacements  in  the  type  II  laminate 
can  be  calculated  from  a  different  set  of  differential  equations  given 
in  reference  [16]  but  are  too  lengthy  to  be  listed  here  along  with  the 
general  solutions. 

It  can  now  be  shown  hov/  moisture  affects  the  crack  spacing  as 
defined  in  this  section. 

3.5  Influence  of  Environmental  Conditioning  on  CDS 

The  influence  of  environmental  conditioning  on  the  CDS  for  type  I 
and  type  II  laminates  is  demonstrated  in  this  section  by  showing  how 
absorbed  moisture  changes  the  predicted  crack  spacing  of  the  CDS  as 
defined  in  section  3.4, 

The  crack  spacing  is  characteristic  of  the  laminate  since  the 
constraining  effect  of  the  stacking  sequence  is  used  together  with  the 
ply  and  laminate  properties  to  calculate  the  distance  between  cracks. 
Although  the  shear  lag  model  can  be  conveniently  explained  in  terms  of 
shear  stress  transfer  of  adjacent  plies  and  the  redistribution  of  stress, 
tile  shear  lag  model  is  more  general  v/hen  formulated  in  terms  of 
normalized  displacement,  u.  Stresses  existing  in  a  cracked  ply  are 
unique  to  the  material  type  whereas  solutions  of  the  shear  lag  model, 
when  calculated  in  terms  of  normalized  displacement,  are  more  general 
and  can  be  applied  to  any  material  type.  For  these  reasons  it  is 
more  general  to  define  the  distance  between  cracks  in  termiS  of  the 
normalized  distance,  x,  when  the  normalized  strain,  du/dx,  reaches  a 
value  of  0.999,  Since  the  crack  spacing  is  more  generally  defined  in 
terms  of  normalized  strain,  du/dx,  and  distance,  x,  it  follows  from 
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the  formulation  in  section  3.4  that  the  effect  of  moisture  absorption 
on  X,  when  du/dx  is  equal  to  0.999,  can  only  be  accounted  for  by  a 
reduction  in  ply  and  laminate  elastic  properties.  A  wet  and  dry 
distance,  x,  at  du/dx  =  0.999  is  calculated  using  wet  and  dry  elastic 
properties,  respectively,  as  listed  in  Table  6.  A  decrease  of  0.4 
percent  is  calculated  for  the  normalized  distance,  x,  between  cracks 
when  1.2  percent  moisture  is  absorbed.  When  the  normalized  distance, 

X,  is  renormalized  by  using  equation  (26),  the  differences  in  wet  and 
dry  distance,  X,  between  cracks,  as  shown  in  Figure  13  is  increased  by 
a  net  9  percent.  The  redistribution  of  stresses  plotted  in  Figure  i3 
is  calculated  by  differentiating  equation  (25)  with  respect  to  X  and 
multiplying  the  resulting  strain  by  E^q  such  that  the  unaisturbed 
stress  in  the  cracked  ply  which  has  recovered  99.9  percent  of  its 
value  is  equated  to  the  ply  strength.  As  a  result  of  these  calculations, 
variations  in  the  ply  strength  due  to  wet  cr  dry  conditioning  have  no 
effect  on  the  predicted  crack  spacing  model  described  in  section  3.4. 

The  net  increase  in  crack  spacing  due  to  moisture  absorption  is  mainly 
due  to  the  factor  shown  in  Equation  (26)  since  the  resin 

rich  shear  transfer  layer  modulus,  is  reduced  more  by  moisture 

than  is  the  laminate  modulus,  E^.  Predicted  crack  spacing  for  wet  and 
dry  type  I  and  type  II  laminates  is  summarized  in  Table  9. 

As  a  final  note,  the  residual  stress  state  in  no  way  affects  the 
rate  of  recovery  of  the  stress  from  the  crack  to  the  undisturbed 
value  of  the  ply  stress  (ply  strength).  The  compressive  wet  residual 
90°  ply  stress,  u^WET’  calculated  in  section  3.2,  results  in  a  higher 
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laminate  load  to  initiate  FPF  even  though  the  transverse  strengtn  is 
reduced  by  moisture  absorption.  After  the  residual  wet  compressive 
stress  has  been  eliminated  by  the  laminate  tension  load  followed  by 
FPF  ♦'he  stress  recovery  from  the  cracked  surface  is  unaffected  by  any 
previous  residual  stress  state.  As  a  result,  the  distance  between 
cracks  is  unaffected  by  residual  stress  states. 

The  predicted  FPF  laminate  loads  and  crack  spacing  discussed  ’n 
this  section  e-d  previous  sections  can  now  be  compared  with  respect  to 
the  damage  at  the  laminate  edge  which  is  recorded  using  the  replica 
tGchnique  discussed  in  section  2.3. 


IV.  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


4.1  Results  of  Environmental  Conditioning 

As  discussed  in  section  2.1  all  specimens  listed  in  Table  1  were 
either  dried  until  all  diffused  moisture  was  desorbed  (Dry)  or 
exposed  to  a  9S%  RH  at  70°C  until  the  maximum  amount  of  diffused 
moisture  was  absorbed  (Wet).  Replicas  of  specimens  designated  with 
asterisks  in  Table  1  were  taken  before  and  after  environmental 
conditioning.  No  damage  was  observed  on  replicas  as  a  result  of 
environmental  conditioning.  Initial  damage  prio?’  to  environmental 
conditioning  was  observed  only  for  specimens  41  (3,4)  as  four  transverse 
cracks  in  the  center  -45"  plies  of  the  type  II  laminates.  When 
specimens  41  (3,4)  v;ere  dried  as  previously  described  no  additional 
damage  resulted  from  the  desorption. 

Before  continuing  with  the  discussion  of  damage  events  in  the 
following  sections,  a  clarification  should  be  made  on  interpreta¬ 
tion  of  damage  as  observed  on  replicas.  The  laminate  edge  damage 
when  transfe'-red  onto  the  acetate  strip  is  assumed  to  be  negligibly 
influenced  by  the  acetone  due  to  the  short  time  of  exposure.  Seven 

succesive  replicas  were  taken  while  the  laminate  load  was  held 

FPF 

constant  at  a  value  greater  than  N^^  and  no  additional  damage  was 
observed  to  have  occurred  due  to  possible  weakening  of  material  at  the 
free  edge  surface  during  the  exposure  to  acetone.  The  image  of  edge 
damage  as  transferred  onto  the  acetate  is  not  as  clearly  seen  ivhen 
compared  to  image  of  edge  damage  when  directly  viewed  through  a 
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microscope.  For  this  reason  the  enlargement  of  a  replica  from  the 
surface  of  fiber  ends  in  a  90°  ply  was  used  in  Figure  3  to  demonstrate 
the  loss  of  resolution.  A  trade  off  is  made  between  quickly  recording 
edge  damage  for  later  observation  as  opposed  to  a  more  time  consuming 
but  more  detailed  microscopic  observation.  Kim  and  Hahn  [11]  success¬ 
fully  used  acoustic  emissions  (AE)  to  record  damage  in  the  form  of  FPF 
but  AE  is  limited  in  this  respect.  The  experimenter  must  ascertain 
whether  the  resolution  of  damage  as  transferred  onto  an  acetate  strip 
has  sufficient  clarity  to  demonstrate  the  damage  of  interest.  In  this 
effort  only  damage  in  the  form  of  cracks,  one  ply  thickness  in  length 
is  of  interest.  Any  cracks  smaller  than  this  are  not  considered  to 
have  eventful  influences  on  the  damage  as  defined  by  the  CDS.  As 
discussed  in  the  introduction,  cracks  most  likely  originate  at  the 
fiber  matrix  interface;  but  until  the  cracks  grow  to  the  size  of  a 
ply  thickness  the  cracks  smaller  than  a  ply  thickness  will  have  no 
effect  on  the  CDS,  as  defined  in  section  3.4. 

4.2  Results  of  Unidirectional  Tests 

The  wet  and  dry  [Og]  and  [90g]  specimens  listed  in  Table  1  were 
loaded  to  failure  as  described  in  section  2.2.  Young's  moduli  for 
both  types  of  specimens,  in  the  wet  and  dry  conditions  are  listed  in 
Figure  15.  Wet  and  dry  strengths  for  both  types  of  specimens  are 
listed  in  Figure  14,  and  compared  graphically  with  other  lamina  and 
laminate  strengths  in  Figure  16.  The  wet  and  dry  v;eakest  link 
strengths  for  the  [90g]  specimens  are  listed  in  Table  10.  The 
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TABLE  10.  SUfWRY  OF  [90]g  WET  AND  DRY  WEAKEST  LINK  STRENGTHS 


Break 

No. 

Dry /Wet 
Average 
Strength  (ksi) 

Wet  -  Dry 
Percent 
Difference 

Dry/Wet 

t  Std.  Dev.  (ksi) 

Dry/ Wet 
No.  of 
Samples 

[Ml 

6.47/4.47 

-27 

0.79/1.02 

4/6 

2 

7.34/4.98 

-32 

.  0.64/0.48 

4/6 

3 

7.65/4.70 

-39 

0.60/0.71 

4/6 

4 

7.92/5.14 

-35 

0.68/0.79 

4/6 

5 

7.68/4.66 

-39 

0.36/0.58 

4/6 

6 

7.39/5.13 

-31 

0.10/0.75 

3/6 

7 

8.17/5.03 

-38 

-  /0.87 

1/6 

8 

7.13/5.46 

-23 

-  /0.06 

1/2 

9 

-  /5.32 

-  /  - 

-/I 

determined  as  0.313  and  0.315 


dry  and  wet  Poisson's  ratios, 
respectively. 

Neither  the  wet  nor  dry  elastic  properties  determined  from  these 
tests  were  used  in  Chapter  3  for  predicting  the  laminate  stress  state. 
Properties  representative  of  experimental  changes  in  properties  due 
to  wet  and  dry  conditions  where  chosen,  as  described  in  section  3.1, 
such  that  the  experimental  wet  value  of  E-j  was  used  as  a  basis  for 
choosing  wet  and  dry  elastic  properties.  These  large  variations 
between  wet  and  dry  values  for  E^,  as  discussed  in  section  3.1,  are 
reproduced  for  the  [Og]  tests  as  shown  in  Figure  15.  A  reduction  of 
9.^  percent  in  wet  transverse  modulus,  E2.  from  the  dry  state,  as 
shown  in  Figure  15,  is  also  representative  of  experimental  changes 

for  E2  as  discussed  in  section  3.1. 

Although  the  experimental  elastic  properties  determined  from  the 
unidirectional  tests  are  not  used  in  the  stress  analysis,  the  wet  and 
dry  strengths  are  used,  as  described  in  section  3.2,  to  predict  the 
laminate  load  required  for  FPF, 

The  vjeakest  link  strength  tests  for  the  dry  [90g]  specimens  appear 
to  reach  an  equilibrium  value,  as  shov;n  in  Figure  14,  which  is 
noticably  higher  than  the  first  dry  failure.  This  same  trend  is  not 
reproduced  in  wet  strengths  which  is  also  graphically  demonstrated  .n 
Figure  14.  The  dry  and  wet  [90gl  strengths  are  graphically  referenced 
with  respect  to  the  laminate  stress  by  using  laminate  theory  as 
shown  in  Figure  16.  It  is  tempting  to  relate  the  experimentally 

16  with  the  first  failure  measured  for 


observed  FPF  plotted  in  Figure 
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the  [90g]  test.  Although  the  graphical  correlation  of  strengtiis  for 
[Og],  [90g],  type  I  and  type  II  tests,  as  shown  in  Figure  16,  is 
instructive,  the  correlation  of  individual  failure  events  should 
not  be  assumed  since  the  [90g]  and  [Og]  tests  are  not  in  the  con¬ 
strained  state  as  are  the  90°  and  0°  plies  of  a  type  I  or  type  II 
laminate.  None  the  less,  using  the  maximum  stress  theory,  as 
described  in  section  3.2,  justifies  comparing  the  [90g]  strengths  and 
experimental  FPF  as  shown  in  Figure  16.  There  are  obvious  differences 
in  the  wet  FPF's  and  first  wet  [90g]  strength  which  obviously  cannot 
be  justified  on  the  basis  of  maximum  stress  theory.  This  would  imply 
that  the  edge  stresses,  which  are  larger  than  the  interior  stresses, 
may  be  responsible  for  the  lov/er  FPF's  shown  in  Figure  16.  These 
differences  will  be  discussed  in  detail  in  section  4.5. 

4.3  Results  of  Static  Tests 

4.3.1  Preliminary  Quasi-Isotropic  Tension  Tests 
As  discussed  in  the  Introduction,  previous  studies  [13,16]  have 
demonstrated  that  damage  which  develops  along  the  edge  of  a  type  I 
laminate  is  different  than  the  damage  which  develops  along  the  edge 
of  a  type  II  laminate.  For  wet  type  I  and  type  II  laminates  it  was 
shown  in  section  3.2  that  the  predicted  residual  stresses  increase  the 
laminate  load  required  to  produce  FPF.  These  differences  in  the  pre¬ 
dicted  wet  and  dry  stresses  existing  in  the  laminate  inte»'ior  or  along 
the  edge  are  significant  as  demonstrated  in  sections  3.2  and  3.3  and 
are  sunnarized  in  Tables  7  and  8.  Corresponding  differences  in  edge 
damage  should  result  from  these  differences  in  wet  and  dry  stresses. 
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The  objectives  of  these  static  tension  tests  are  to  initially 

isolate  unique  free  edge  damage  states  in  wet  and  dry  type  I  and  type 

II  laminates  and  show  how  these  unique  damage  states  develop  into  the 

CDS  when  the  laminate  load  is  increased.  Specimens  listed  in  Table  1 

were  tested  as  described  in  section  2.3.  The  results  of  these  four 

tension  tests  are  listed  in  Tables  9  and  11. 

As  shown  in  Table  11  there  is  a  substantial  difference  between 

the  predicted  and  observed  damage  in  the  90°  plies.  The  predicted 

FPF 

values  of  a  ,  a  and  Ny  v;hich  v;ere  calculated  in  sections  3.2  and 

X  Z  A 

3.3  assumed  a  damage  free  laminate  with  all  plies  having  the  same 

thickness.  Only  trends  in  damage  events  are  predicted  in  Table  11 

FPF 

without  reference  to  predicted  values  of  a^,  and  Ny  since  the 
90°  plies  in  the  type  I  laminate  vjere  only  one  ply  thickness  and  a 
large  amount  of  initial  damage  in  the  type  II  laminates  existed  in 
the  form  of  longitudinal  and  transverse  cracks.  For  instance,  the 
trend  in  the  damage  predicted  for  type  I  wet  specimen,  SI 2,  assumed 
delaminations  v/ould  occur  as  the  first  damage  event  followed  by 
transverse  cracking  since  is  much  larger  than  at  the  edge 
(Y/B  =1.0).  As  a  result  of  these  irregularities  a  detailed 
comparison  betv/een  predicted  and  observed  damage  is  not  discussed 
in  this  section.  The  necessary  comparisons  between  predicted  and 
observed  damage  is  discussed  in  the  following  section  where  the 
test  specimens  were  fabricated  v/ith  the  correct  ply  thickness  and 
no  initial  damage. 

Although  the  differences  in  the  predicted  damage  and  experimental 


TABLE  n.  DAMAGE  DEVELOPMENT  INTO  CHARACTERISTIC  DAMAGE  STATE  FOR 
PRELIMINARY  QUASI-ISOTROPIC  TENSION  TESTS 
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observations  are  most  likely  due  to  the  iri-egularities  in  fabrication, 
the  damage  initially  isolated  in  each  test  is  unique.  The  growth  of 
these  initially  unique  damage  states  resulted  in  the  same  CDS  in  the 
laminate  whether  wet  or  dry.  The  final  CDS  prior  to  failure  is 
suimarized  in  Table  11,  where  the  load  level,  (LL),  and  description 
of  damage  near  the  fractured  region  are  listed  in  an  abbreviated 
format.  Except  for  specimen  S22,  all  damage  leading  to  the  final 
fracture  was  observed  to  grow  from  newly  formed  damage  while  under 
load.  The  final  fracture  surface  of  the  S22  specimens  was 
observed  to  grow  from  initial  damage  which  was  observed  along  the 
free  edge  in  the  90°  and  45°  plies  prior  to  the  load.  The  damage, 
whether  induced  by  the  load  or  existing  prior  to  the  load,  does  not 
change  the  crack  pattern  of  the  CDS  prior  to  fracture. 

Contrary  to  the  irregularities  (ply  thickness,  initial  damage) 
discussed,  the  preliminary  tests  provided  evidence  in  the  early  part 
of  the  experimental  program  that  unique  damage  states  followed  by 
growth  into  the  CDS  with  increasing  load  could  be  observed  by  using 
the  replica  technique. 

While  under  a  quasi-static  tension  load  all  damage  eventually 
grew  into  the  same  crack  patterns  which  were  characteristic  of  that 
laminate  (CDS),  '’’hese  patterns  grew  independently  of  the  uniqueness  of 
the  first  formation  of  damage.  Even  if  the  growth  originated  from 
initial  damage  due  to  fabrication  there  was  no  appreciable  difference 
in  the  crack  patterns  prior  to  failure.  From  the  results  of  these 
preliminary  tests  the  following  test  and  improvements  were  established; 


Using  the  interactive  coaxing  technique,  described  in  section  2.5, 
attempt  to  coax  the  first  formation  of  damage  (i.e.  initial  damage 
due  to  fabrication)  into  a  final  damage  state  different  than  the  CDS. 

If  the  order  of  occurrence  of  individual  damage  events  (FPF, 
equilibrium  crack  spacing,  edge  delamination)  due  to  quasi-static 
or  cyclic  loads  are  to  be  recorded  then  more  replicas  should  be 
taken  at  smaller  load  intervals. 

4.3.2  Development  of  CDS  Due  to  a  Quasi-Static  Tension  Load 

The  preliminary  tests  demonstrated  that  unique  damage  states  and 
subsequent  gro.vth  could  be  recorded  using  the  replica  technique. 
Unfortunately  the  stress  state  predicted  in  section  3.2  and  3.3  could 
not  be  compared  with  the  damage  events  observed  in  the  preliminary 
tests  due  to  irregularities  detected  in  the  test  specimens  after  fab¬ 
rication.  In  this  section  specimens  with  no  initial  damage  and  correct 
ply  thicknesses  are  tested  as  described  in  section  2.4. 

No  replicas  were  taken  for  the  specimens  listed  in  Table  3  which 
were  quasi-statically  loaded  to  fracture.  The  difference  between 
type  I  and  type  II  laminate  strengths  is  graphically  demonstrated 
in  Figure  14  for  both  wet  and  dry  conditions  which  are  compared  with 
the  experimental  dry  strengths  from  reference  [14].  When  these  vjet 
and  dry  strengths  are  graphically  referenced  with  respect  to  the 
laminate  stress  state,  as  uhown  in  Figure  16,  the  absorbtion  of 
moisture  is  observed  to  have  reduced  the  differences  between  the 
type  I  and  type  II  dry  strengths.  It  is  interesting  to  note  that 
a  similar  reduction  in  scatter  is  observed  for  the  wet  [Og]  tests  also 


shov/n  in  Figure  16.  Since  the  type  I  and  type  li  strengths  are 
dominated  by  the  0°  plies  these  similarities  in  wet  and  dry  strengths 
are  not  surprising.  Moduli  determined  for  the  wet  and  dry  type  I 
and  type  II  laminates  are  listed  in  Figure  15. 

As  discussed  in  the  introduction  the  experimental  differences 
between  type  I  and  type  II  laminate  strengths  provide  a  basis  for 
investigating  how  the  damage  develops  under  load  and  influences  the 
laminate  strength.  As  previously  noted  in  section  3.4,  unoerstanding 
how  the  different  CDS  crack  patterns  influence  the  distribution  of 
stresses  prior  to  fracture  is  a  formidable  task  which  is  still  under 
investigation.  At  best  the  effort  in  this  section  is  t.o  demonstrate 
how  the  unique  damage  states  due  to  the  wet  conditioning  a'^fect  the 
CDS  as  defined  in  section  3.4. 

The  last  four  specimens  listed  in  Table  3  were  tested  as  described 
in  section  2.4  such  that  first  formation  of  damage  and  the  subsequent 
growth  can  be  studied  and  compared  with  predictions  for  FPF  and  crack 
spacing  from  section  3.2  and  3.4  respectively.  The  results  of  these 
tests  are  briefly  summarized  in  Tables  9  and  12. 

The  predicted  value  for  and  stresses  in  the  90°  ply  as 

listed  in  Table  12  are  taken  from  Table  7.  The  predicted  from 

section  3.2,  are  also  listed  in  Table  12  as  the  laminate  load,  (LL), 
shown  next  to  the  T  designation  for  90°  ply  transverse  crack. 

Although  no  predictions  for  the  laminate  load  are  made  for  edge 
delaminations,  the  predicted  and  stress  would  indicate  that 

no  del aminotions would  occur  at  a  lower  laminate  load  for  the  wet 
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conditioned  type  I  laminate. 

Although  no  elevated  temperature  tests  were  included  in  the 
experimental  program,  the  trend  at  higher  test  temperatures  for  type  I 
wet  laminates  is  to  cause  the  predicted  value  of  to  decrease  more 
than  the  value  predicted  for  as  temperature  is  increased  with  the 
laminate  load  held  constant.  This  trend  would  obviously  increase  the 
laminate  load  required  to  cause  FPF  but  would  lower  the  laminate  load 
required  to  cause  delaminations.  The  effect  of  increasing  test 
temperatures  is  similar  to  the  effect  of  swelling  due  to  moisture 
absorbtion.  Higher  test  temperatures  stress  relieve  the  curing 
stresses.  This  results  in  a  larger  predicted  residual  value  for  the 
compressive  wet  90°  ply  stress,  It  may  be  possible  to  verify 
this  trend  experimentally  by  testing  type  I  wet  laminates  at  elevated 
temperatures  and  show  edge  delaminations  occurring  before  transverse 
cracking.  Although  no  elevated  temperature  tests  were  performed  the 
swelling  of  the  type  I  laminate  due  to  moisture  experimentally 
verifies  this  trend.  The  observed  first  occurrence  of  a  wet  transverse 
crack  (FPF)  occurs  between  1100  Ib/in  and  1200  Ib/in  which  is  con¬ 
siderably  higher  than  the  first  transverse  cracks  observed  in  the 
interval  (500-600)  Ib/in  for  the  dry  case.  Also  the  laminate  load 
required  to  cause  delamination  has  dropped  from  (1800-2000)  Ib/in  for 
the  dry  case  to  (1100-1200)  Ib/in  for  the  wet  case.  Interestingly 
enough  both  wet  transverse  cracks  and  wet  delaminations  appear  to 
occur  within  the  same  replica  load  interval  of  (1100-1200)  Ib/in,  as 
shown  in  Figure  17. 
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A  similar  trend  for  the  type  II  wet  and  dry  laminates  is  not 
experimentally  demonstrated  as  expected  since  the  predicted  value  of 
2.87  ksi  for  in  a  dry  type  II  90°  ply  is  larger  than  the  wet 
value,  1.71  ksi.  Both  v/et  and  dry  values  for  are  much  lower  than 
at  FPF  and  therefore  delaminations  are  not  likely  to  occur  for  the 
type  II  wet  cr  dry  laminates,  as  is  experimentally  verified  in  the 
replicas  shown  in  Figure  18  and  19  which  were  taken  prior  to 
fracture. 

The  equilibrium  crack  spacing  in  the  wet  90°  plies  for  both  type 
1  and  type  II  laminates  is  not  affected  by  the  compressive  wet 
residual  stress,  a^.  As  discussed  in  section  3.5  the  residual 
stresses  in  the  wet  laminates  only  result  in  increasing  the  laminate 
load  required  for  FPF  and  do  not  alter  the  crack  spacing  of  the  CDS 
as  defined  in  section  3.4.  The  results  of  the  crack  spacing  are 
compared  with  the  predicted  values  listed  in  Table  9.  Due  to  the 
simultaneous  occurence  of  delaminations  and  transverse  cracks. in  the 
type  I  wet  90°  plies  the  transverse  cracks  are  not  as  clearly  defined 
as  in  the  dry  condition.  The  wet  type  I  90°  ply  crack  shown  in 
Figure  17  appears  to  be  less  vertically  inclined  than  the  dry  transverse 
cracks  and  tend  to  interact  with  the  delaminations  resulting  in  an 
equilibrium  crack  spacing  which  is  twice  as  large  as  the  dry 
equilibrium  crack  spacing  as  shown  in  Table  9.  Both  wet  and  dry 
equilibrium  crack  spacing  are  listed  in  Table  9  and  are  compared  with 
predicted  crack  spacing  from  section  3.5.  In  summary  the  change  in 
the  order  of  damage  events  due  to  moisture  can  obviously  affect  the 
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FIG’JRF.  19.  FUl.l.Y  DLVELOPfD  CHARACTER  I  STIC  (JAHAGE  STATE  FOR  TYPE 
II  WET  LAMINATE  WITH  N  =34CO  LB/IN 


equilibrium  crack  pattern  in  the  90°  p'y  for  the  type  I  laminate. 

Although  the  simultaneous  occurrence  of  delaminations  and  transverse 
cracks  is  the  cause  of  the  large  crack  spacing,  the  final  pattern  of 
cracks  recorded  at  2200  Ib/in  for  the  type  I  dry  and  at  2400  Ib/in  for 
the  type  I  wet  are  strikingly  similar  as  shown  in  Figures  20  and  21, 
respectively.  Small  hair  like  strands  extend  from  the  acetate  surface 
when  the  replica  is  removed  before  the  acetone  has  evaporated.  When 
the  acetate  strip  is  sandwiched  bet'.veen  glass  plates  for  viewing, 
these  hair  like  strands  are  seen,  as  shovjn  in  Figure  21,  as  cracks 
extending  from  the  delaminaticn  cracks.  Even  though  the  90°  ply  crack 
pattern  is  substantially  altered  in  the  type  I  laminate  due  to 
moisture  absorption,  the  final  crack  pattern  prior  to  fracture  is 
not  significantly  altered.  If  the  characteristic  wet  and  dry  crack 
patterns  are  similar,  it  follows  that  the  influence  of  the  wet  and 
dry  damaged  stress  state  on  the  laminate  strength  should  also  be 
similar.  This  observation  is  experimentally  verified  in  Table  14 
where  the  average  type  I  laminate  dry  strength  is  increased  by  only 
4  percent  when  1.2  percent  moisture  is  absorbed.  Similarly,  the 
type  II  wet  and  dry  crack  patterns  prior  to  fracture  are  not  sub¬ 
stantially  different  as  shown  in  Figures  18  and  19.  In  summary,  the 
effect  of  moisture  can  significantly  alter  crack  patterns  in  the 
type  I,  90°  plies  but  the  effect  on  the  final  CDS  and  subsequent 
laminate  strength  is  negligible. 

The  cracks  in  the  45°  plies  were,  not  included  when  predicting 
the  crack  pattern  of  the  CDS  as  defined  in  section  3.4.  The  45°  ply 
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FIGURK  20.  FULLY  DEVELOPED  CHARACTERISTIC  DAMAGE  STATE  FOR  TYPE  I 
DRY  LAMINATE  N  =2200  LB/IN 
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crack  spacing  predicted  by  the  shear  lag  model  [16]  assumes  that  the 
45°  ply  crack  patterns  develop  independently  of  the  90°  ply  cracks. 

For  wet  or  dry  type  I  and  type  II  laminates,  nearly  all  45°  ply 
cracks  which  develop  at  the  laminate  edge  appeared  to  be  influenced 
by  transverse  cracks  in  the  adjacent  90°  plies,  as  shown  in  Figures 
18  through  21 . 

4.4  Results  of  Fatigue  Tests 

4.4.1  Coaxing  Out  Different  CDS 

Four  specimens  listed  in  Table  4  were  cyclically  loaded  using  the 
interactive  coaxing  technique  described  in  section  2.5.  The  basis  for 
this  series  of  tests  v^as  suggested  after  observing  the  growth  of 
initial  damage  along  the  edge  of  the  preliminary  quasi-isotropic 
tension  tests  of  section  4.3.1. 

As  previously  discussed  in  section  2.5,  there  is  no  clearly 
defined  experimental  procedure  for  the  coaxing  fatigue  tests.  The 
material  presented  in  Chapter  III  which  has  been  compared  with 
experimental  observation  in  sections  4,3.1  and  4.3.2  provides  a  better 
basis  at  this  time  on  wiiich  to  outline  the  following  guidelines.  As 
described  in  section  4.3.1,  regardless  of  the  uniqueness  of  the  first 
formation  of  damage,  all  damage  developed  into  the  same  CDS  prior  to 
fracture.  This  point  was  reinforced  in  section  4.3.2  as  shown  in 
Figures  18  through  21.  The  objective  of  the  coaxing  tests  is  not  to 
observe  the  normal  sequence  of  damage  events  but  to  preload  the 
specimens  until  a  unique  damage  state  is  observed  followed  by  a  cyclic 
load  with  a  maximum  load  equivalent  to  or  less  then  the  preload.  The 


microfiche  card  reader  was  used  to  allow  for  interpretation  of  the 
damage  during  the  test  preload.  The  cyclic  load  which  would  peak  at 
or  below  the  prelead  level  would  give  the  existing  unique  damage  an 
opportunity  to  grow  into  a  damage  state  different  than  the  normal 
sequence  of  damage  events  demonstrated  in  section  4.3.2. 

Each  of  the  four  laminates  listed  in  Table  4  were  preloaded  as 

described  in  Table  13.  The  guidelines  for  each  of  the  four  tests  are 

also  listed  in  Table  13.  Except  for  the  type  I  wet  test,  the  guide¬ 
lines  for  coaxing  the  growth  of  new  damage  are  reasonable  v^hen 
considering  the  normal  sequence  of  damage  events  as  d’seussed  in 
section  4.3.2.  The  inherent  nature  of  the  large  equilibrium  crack 
spacing  for  the  type  I  wet  laminate  was  not  understood  when  the  coaxing 
fatigue  tests  were  being  conducted.  As  a  result,  the  higher  2000  Ib/in 
laminate  preload  was  applied  and  subsequently  cycled.  A  more 
reasonable  preload  of  1100  Ib/in  with  a  maximum  cyclic  load  of  1000 

lbs  would  have  been  a  wiser  choice.  This  lov/er  cyclic  load  may  have 

allowed  the  initial  delaminations  to  grow  independently  of  the 
initial  transverse  cracks.  As  a  result,  the  lower  cyclic  load  may 
have  produced  an  even  larger  equilibrium  crack  spacing  than  recorded 
in  Table  9  for  the  type  I  wet  laminate. 

The  type  II  wet  laminate  was  preloaded  and  cycled  at  a  load  such 
that  no  damage  could  be  observed  after  the  preload.  Although  no 
initial  damage  due  to  fabrication  was  pr'esent,  the  growth  of  micro¬ 
scopic  damage  in  the  form  of  randomly  spaced  transverse  cracks  was 
idealized;  but,  as  shown  in  Table  13,  this  unique  CDS  never 


TABLE  13.  DAMAGE  DEVELOPMENT  INTO  CMARACTERISTIC  DAMAGE  STATE 
FOR  THE  COAXING  FATIGUE  TESTS 
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matprialized,  even  after  2.5  minion  cycles.  The  only  coaxing 
fatigue  test  that  demonstrated  a  damage  event  different  from  the 
normal  sequence  was  the  type  II  dry.  The  type  II  dry  specimen  was 
preloaded  such  that  equilibrium  spacing  was  realized  with  no 
delaminations,  ND,  and  no  cracks  in  the  45°  ply.  No  guideline  was 
set  for  the  cyclic  portion  of  the  test.  More  45°  ply  cracks,  as 
shown  in  Figui-e  22,  appeared  to  grow  independently  of  the  90  cracks 

in  the  adjacent  ply. 

Regardless  of  the  guidelines  used  for  coaxing,  no  large 
differences  in  the  crack  patterns  prior  to  fracture  were  observed. 
Except  for  the  type  I  wet  test  the  residual  strengths  as  listed  in 
Table  13  are  not  significantly  different  from  the  residual  strengths 
listed  in  Figure  16. 

4.4.2  Development  of  CDS  Due  to  a  Cyclic  Load 

The  coaxing  cyclic  load  described  in  the  previous  section 
demonstrated  an  attempt  to  coax  the  growth  of  damage  which  was  dif¬ 
ferent  than  the  normal  sequence  observed  in  the  static  tension  tests 
In  this  section  the  normal  sequence  of  damage  events  due  to  a  steady 

state  cyclic  load  is  demonstrated. 

All  steady  state  fatigue  tests  listed  in  Table  4  were  tested  as 
described  in  section  2.5.  All  replica  specimens  were  preloadec  to 
2000  lb  in  200  lb  load  increments  followed  by  a  10  Hz  cyclic  load 
at  R  =  0.1.  All  remaining  steady  state  fatigue  tests  which  were  nov 
replicated  were  preloaded  to  50  ksi  followed  by  a  10  Hz  cyclic  load 
at  R  =  0.1.  Replicas  and  moduli  were  recorded  during  the  steady 
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NOTE:  All  other  symbols  defined  in  Table  11. 


state  fatigue  tests  as  outlined  in  Table  4.  Results  of  the  replica 
steady  state  fatigue  tests  are  summarized  in  Table  14.  Results  of 
the  steady  state  fatigue  tests  with  no  replicas  are  summarized  in 
Figures  14  and  16. 

A  brief  summary  of  the  edge  damage  development  as  observed  on 
replicas  due  to  the  preload  followed  by  a  steady  state  cyclic  load  is 
presented  in  Table  14.  During  the  preload  of  type  I  dry  and  wet 
tests,  transverse  cracking  and  delaminations  occurred  at  laminate 
load  levels  (LL)  similar  to  LL  reported  for  the  static  tension  tests 
in  Table  12.  Only  2000  cycles  (2K)  were  required  to  grow  the 
initially  formed  delaminations  into  a  single  delamination,  "full", 
running  the  length  of  the  laminate.  Although  the  higher  predicted 
value  for  the  wet  as  discussed  in  section  4.3.2  is  responsible  for 
delaminations  occurring  sooner  during  the  preload,  these  same 
arguements  can  not  be  used  to  predict  that  the  v/et  delaminations 
should  grow  faster  when  cyclically  loaded.  For  both  wet  and  dry  type 
I  laminates,  delaminations  grow  into  the  fully  developed  state  after 
the  same  number  of  cycles.  This  would  imply  that  the  damaged  stress 
state  differs  considerably  from  the  undamage  state  with  reference 
to  near  the  edge. 

The  damage  resulting  from  preloading  wet  and  dry  type  II 
laminates  occurs  at  loads  similar  to  those  reported  in  Table  12.  No 
delaminations  occur  in  the  type  II  90°  plies  during  the  preload  or 
cycling.  Transverse  cracking  in  wet  or  dry  45°  plies  of  the  type  II 
laminates  results  from  the  cyclic  load  and,  in  most  cases  observed, 
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occurs  after  the  same  number  of  cycles.  In  general,  delaminations 
betvveen  the  +45  and  -45°  plies  of  the  type  H  wet  laminate  occurs 
after  a  larger  number  of  cycles  when  compared  to  the  type  11  dry 
laminate.  Damage  in  the  wet  and  dry  type  I  45°  plies  was  not 
recorded  due  to  largo  delaminations  which  occurred  at  approximately 
2500  lbs.  These  large  delaminations  did  not  allow  for  accurate 
replication  of  edge  damage. 

The  final  and  most  important  observation  was  that  the  final 
recorded  damaged  crack  patterns  of  the  CDS  were  not  significantly 
altered  by  the  wet  or  dry  conditioning.  As  a  result  the  residual 
strengths  of  the  nonrepl icated  test  specimens  listed  in  Figure  16 
demonstrate  negligible  variations  due  to  wet  or  dry  conditioning. 

A  summary,  of  results  for  the  nonrepl icated  steady  state  fatigue 
tests  are  listed  in  Tables  14  and  15.  Not  only  are  the  differences 
small  between  wet  and  dry  residual  strength  after  100  K  cycles  but, 
as  compared  in  Figure  14,  the  difference  in  residual  and  static 
strengths  are  small.  These  same  observations  are  not  true  for  changes 
in  static  or  cyclic  moduli  due  to  wet  or  dry  conditioning,  as  shown 
in  Figure  15.  The  trend  after  a  100  K  cycles  is  demonstrated  as  a 
larger  change  of  moduli  for  the  wet  and  dry  type  I  laminates,  with 
the  smallest  change  in  moduli  occurring  after  100  K  cycles  for  the 
type  II  dry  laminate. 

4.5  Comparison  of  Stress  Analysis  with  Experimental  Observations 

As  shown  in  Table  12  and  Figure  16  there  are  obvious  differences 
between  dry  and  wet  experimental  laminate  loads  (N  )  required  to 
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produce  the  first  transverse  cracks  in  the  90°  plies  as  observed  on 
the  acetate  strips.  Although  the  data  is  scattered  due  to  the  load 
increments  over  which  replicas  are  taken,  a  trend  is  evident  between 
the  wet  and  dry  type  I  and  type  II  laminate  loads  (N^)  required  to 
produce  First  Ply  Failure  (FPF)  in  the  90°  plies.  The  difference 
betv^een  ^or  type  I  wet  and  type  II  wet  is  greater  than  the 
difference  between  for  type  I  dry  and  type  II  dry.  Investigations 
by  Kim  and  Hahn  [11]  have  shown  differences  between  the  dry  and  wet 
type  I  laminate  FPF  load  but  these  were  attributed  only  to  the 
interior  residual  stress  state  as  predicted  by  two-dimensional 
laminate  plate  theory. 

It  is  particularly  interesting  that  the  differences  between  type  I 
and  type  II  dry  laminate  FPF  loads  is  larger  for  the  same  laminates  in 

the  wet  condition.  This  trend  is  not  surprising  when  we  recall  the 

through  the  thickness  and  stress  distributions  near  the  laminate 
free  edge  (Y/B=0.998)  as  shown  in  Figure  12.  In  Figure  12  there  is  an 

obvious  difference  between  the  interior  and  edge  and  stresses  in 

the  90°  plies  for  the  four  laminates.  The  differences  between  the 
interior  and  edge  and  0^  stresses  for  the  type  I  wet  laminate  is 
much  greater  than  the  difference  between  the  same  interior  and  edge 
stresses  in  the  type  II  wet  which  is  at  the  same  laminate  load 
The  difference  in  90°  ply  interior  and  edge  stresses  for  type  I  wet 
laminate  is  also  larger  than  the  differences  shown  for  the  type  I 
and  type  II  dry  laminates  which  are  at  lower  laminate  loads. 

These  predicted  trends  in  the  stress  concentrations  of  0^  and  0^ 
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near  the  laminate  edge  as  shown  in  Figure  1?  could  be  used  as  a  basis 

for  explaining  the  differences  previously  noted  in  the  experimental 

data.  It  is  tempting  to  use  these  predicted  edge  stress  states  together 

with  some  failure  criterion  in  order  to  explain  the  experimental 
FPF 

differences  in  as  previously  discussed.  Values  of  the  free 
edge  stresses  predicted  by  finite  element  models  are  approximations. 
Crossman  [10]  points  out  that  only  trends  in  the  predicted  stress 
state  can  be  utilized  since  approximations  are  made  for  properties 
in  the  lamina  2-3  plane.  Also  the  magnitude  of  is  questionable 

at  the  free  edge,  Y/B=1.0,  since  the  composite  material  can  no  longer 
be  modeled  as  behaving  as  a  homogeneous  material  one  or  two  fiber 
diameters  from  the  free  edge.  With  these  peculiarities  in  the  stress 
state  established,  it  is  with  great  caution  that  the  author  proposes 
the  following  model  which  will  be  used  to  predict  FPF  in  90®  plies  for 
any  laminate. 

The  proposed  model  assumes  that  first  ply  failure  initiates  in  the 
■90°  ply  near  the  laminate  free  edge  since  the  values  predicted  for 
and  in  the  90°  ply  are  larger  near  the  free  edge  then  in  the 
laminate  interior.  Although  more  accurate  predictions  of  and 
stresses  are  obtained  in  this  dissertation  by  using  better  approxima¬ 
tions  for  the  elastic  properties  in  the  lamina  2-3  plane,  these 
improvements  in  predicting  and  stresses  are  still  incomplete 
from  a  continuum  viewpoint  near  the  edge.  A  model  is  needed  which 
can  predict  failure  without  reference  to  an  exact  prediction  of  an 
experimental  stress  near  a  stress  concentration.  Whitney  and  Nuismer 


[41]  point  out  that  determining  the  strength  of  a  material  from  the 
maximum  stress  at  a  point  is  questionable,  especially  when  the  maximum 
stress  is  highly  localized.  Whitney  and  Nuismer  further  show  that  the 
localized  nature  of  a  stress  concentration  near  a  hole  of  radius,  R, 
in  any  quasi-isotropic  laminate  can  be  empirically  modeled  by 
evaluating  a  characteristic  distance,  dQ,  from  the  elasticity 
solution  (34),  where  the  stress  ratio  o^/a  in  equation  (34)  is 
redefined  in  terms  of  experimental  notched,  and  unnotched, 

CTq,  laminate  strengths. 

a^(C)/CT  =  2/(2+C^+3c^)  =  0^,/oq  (34) 

where 

C  =  R(R+dg) 

R  =  Hole  radius 

Oy  =  Stress  near  the  stress  Concentration 
a  =  Far-field  stress 
=  Notched  Strength 
Oq  =  Unnotched  Strength 

dg  =  Characteristic  distance  from  stress  concentration 
In  summary  Whitney  and  Nusimer  demonstrate  that  although  a  critical 
stress  at  some  distance,  dp,  is  the  cause  of  some  lower  laminate 
strength,  the  changes  in  experimental  strengths,  can  be  predicted 

for  different  hole  radii,  R,  without  predicting  an  experimental  value 
for  this  critical  stress. 

Although  these  stress  concentrations  were  modeled  near  holes  in 


quasi-isotropic  laminates,  the  same  philosophy  can  be  extendeu  to 
stress  concentrations  in  individual  plies  near  the  laminate  free  edge. 
The  objective  of  the  proposed  model  is  to  predict  failure  in  the  90 
ply  without  reference  to  predicting  an  experimental  stress  near  the 
stress  concentration.  In  the  90“  ply,  the  stress  concentration  near 
the  laminate  edge  is  not  as  easily  modeled  as  the  stress  concentration 
near  a  hole  in  an  isotropic  plate.  Consequently  there  is  no  closed 
form  solution  predicting  and  in  the  90°  ply  as  a  function  of  the 
distance  from  the  stress  concentration  since  and  are  predicted 
from  a  finite  element  model.  We  do  not  require  that  the  predicted 
stress  state  near  the  laminate  edge  at  some  distance,  Zq,  must  be 
numerically  equal  to  the  actual  stress  state  at  that  point.  We  only 
require  that  the  variations  of  and  at  some  point  near  the  edge 
due  to  variations  in  stacking  sequence  together  with  wet  or  dry 
conditioning  be  realistically  predicted  by  the  finite  element  model. 

The  point  at  which  variations  in  and  are  to  be  predicted  is  not 
arbitrarily  chosen.  As  suggested  by  Crossman  [10]  from  a  continuum 
viewpoint  the  stresses  evaluated  at  five  fiber  diameters  from  the 
edge  is  more  realistic  than  stresses  predicted  7  cm  from  the  edge. 
Although  there  is  no  closed  form  equation  for  evaluating  this  distance, 
the  author  chooses  7.7  fiber  diameters  or  Y/B=0.993  as  a  reasonable 
distance  from  the  edge  to  evaluate  the  FPF  stress  state  in  the  90 
ply  for  wet  or  dry  type  I  and  type  II  laminates.  When  these  predicted 
90°  ply  FPF  stress  states  for  each  case  evaluated  in  this  investigation 
are  substituted  into  the  tensor  polynomial  <^ailure  criterion  polynomial 
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of  Appendix  A  together  with  F^.  and  listed  in  Table  15,  then 
the  surviving  terms  of  equation  (A. 2)  are  written  in  equation  (35). 

The  terms  which  remain  and  contribute  significantly  to  the  failure 
function  (FF)  are  in  agreement  with  the  work  of  Herakovich,  Nagarkar 
and  O'Brien  [42]. 

2ax/T  -  -  (cx/T)^  -  (o^/T)2  =  ff  (35) 

where  the  transverse  strengths  of  the  90®  ply  in  the  lamina  2-3  plane 

are  assumed  equal  *^22~’^33"^'^^^^‘ 

The  model  uses  equation  (35)  to  empirically  evaluate  a  parameter 

T  which  is  representative  of  the  failure  state  for  the  90°  ply  in  any 

wet  or  dry  stacking  sequence.  The  and  stresses  are  predicted  by 

the  FEM  at  Y/B=0. 998  using  the  experimentally  measured  and  the 

failure  function,  ff,  is  assumed  1.0  since  cracks  in  the  90°  ply  are 

experimentally  observed  on  replicas  when  the  experimental  laminate  load 
FPF 

Njj  is  applied  to  the  laminate. 

The  proposed  model  does  not  imply  that  the  and  stresses 
predicted  at  Y/B=0.998  by  the  FEM,  when  an  experimental  is 
applied,  are  the  same  stresses  which  exist  in  the  -.ttaterial  at  Y/B=0.998 
when  FPF  occurs.  Since  the  model  requires  that  the  FEM  accurately 
calculate  variations  in  and  due  to  differences  in  stacking 
sequence  together  with  wet  or  dry  conditioning  then  the  transverse 
strength,  T,  calculated  from  a  and  a_  is  an  effective  parameter,  T 
and  the  proposed  model  is  used  to  show  variations  in  due  to 

A 

variations  in  stacking  sequence  together  with  wet  or  dry  conditioning. 


TABLE  15.  EFFECTIVE  TRANSVERSE  EDGE  STRENGTHS,  LAMINA  STRENGTHS 
AND  DEFTNITIONS  OF  F.  AND  F.. 
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If  the  values  of  and  stresses  predicted  by  FEM  are  the  actual 
stresses  in  the  material  and  Y/B=0.998  is  the  point  where  failure 
occurs,  then  is  indeed  a  mechanical  property  of  the  laminate. 

The  model  used  to  study  trends  and  differences  in  this  dissertation 
ascribes  no  physical  significance  to  but  requires  that  it  be 

a  constant  much  the  same  as  dg  was  empirically  evaluated  as  a 
constant.  With  the  character  of  the  parameter  established  the 
transverse  strength,  T,  in  equation  (35)  becomes  an  effective  parameter 
to  be  solved  for  as  shown  below  in  equation  (36). 

Vf  =  -  °x  ■  4 

where  and  are  stresses  calculated  at  a  dis'^ance  7.7  fiber  diameters 
from  the  laminate  edge  using  the  FEM  described  in  Section  3.3. 

As  noted,  the  proposed  model  requires  that  only  those  terms  listed 
in  equation  (35)  are  numerically  significant  when  calculating  a  value 
for  ff.  Terms  listed  in  (A. 2)  containing  which  are  negligible  near 
the  edge,  become  significant  as  the  edge  stresses  recover  to  the  stress 
state  in  the  interior  of  the  laminate.  Therefore  the  model  as  defined 
in  equation  (35),  is  valid  over  a  small  distance  from  the  edge.  This 
distance,  ag,  cannot  be  smaller  than  is  realistic  from  a  continuum 
point  of  view  and  cannot  be  larger  than  a  certain  value  for  whici.  terms 
other  than  those  listed  in  equation  (35)  significantly  contribute  to  ff. 

where  li  is  the  minimum  distance  over  which  the  composite  material  can 
be  modeled  as  a  continuum  and  1^  is  the  distance  over  which  the  four 
terms  in  equation  (35)  contribute  99.9  percent  of  the  failure  function. 


Obviously  the  lower  limit  is  arbitrarily  based  on  individual 

interpretation  of  the  finite  element  grid  together  with  the  heterogeneous 

nature  of  the  fiber/matrix  material  system.  The  upper  limit  depends  on 

values  assigned  to  the  tensor  polynomial  coefficients  and 

(i»j=l>2,3)  together  with  the  stresses  calculated  from  the  FEM.  For 

the  laminates  studied  in  this  investigation  a  lower  limit  of  five  fiber 

diameters  is  chosen.  In  general  the  upper  limit  for  the  distance,  aQ 

from  the  edge  must  be  less  than  one  laminate  thickness  since  the  edge 

stress  state  has  recovered  to  the  interior  solution  at  this  distance 

and  other  stress  components  contribute  to  FF.  At  Y/B=0.998  or  7.7 

graphite  fiber  diameters  from  the  edge,  the  stresses  calculated  from 

the  FEM  described  in  section  3.3  are  combined  with  coefficients  F.  and 

F..  as  listed  in  Table  15.  For  all  wet  and  dry  laminate  stress  states 
>3 

evaluated  at  Y/B=0.998  all  terms  other  than  those  listed  in  equation  (35) 

contributed  less  than  +0.1%  to  the  final  value  of  ff.  AHhough 

could  be  more  rigorously  defined,  the  value  aQ=l-Y/B=0.002  satisfies 

the  present  definition  expressed  in  the  inequality  (37). 

The  effective  transverse  strengths  listed  in  Table  15  are  calcu- 

lated  as  a  function  of  the  experimental  ‘  laminate  loads  listed  in 

FPF 

Table  12.  The  general  linear  relationships  for  the  laminate  loads 
and  the  transverse  strengths  for  wet  and  dry  type  I  and  type  II 

r  pr 

laminates  are  plotted  in  Figure  23  along  with  the  experimental 
in  order  to  demonstrate  the  trend  discussed  at  the  beginning  of  this 
section.  Although  the  values  calculated  for  T^^^  in  Table  15  are 
scattered  due  to  the  replica  load  increment.  Figure  23  demonstrates 
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that  can  be  used  as  a  parameter  to  predict  changes  in  the 
due  to  stacking  sequence,  where  for  the  wet  case  has  a  lower 
value  than  the  ^eff  for  the  dry  case.  The  model  verifies  that  the 
difference  between  experimental  for  type  I  wet  and  type  II  v/et 
laminates  is  larger  than  the  difference  between  experimental  for 
type  I  dry  and  type  II  dry  laminates.  Therefore  the  trend  is  reproduced 
by  the  model  and  the  mechanism  responsible  for  the  differences  in 
is  due  to  the  variations  in  values  predicted  for  a  and  o  at 

X  2 

Y/B=0.998. 

Figure  24  shows  a  final  comparison  of  the  effect  that  edge  stresses 

FPF 

have  on  the  laminate  load  )  required  to  initiate  cracking  in  the 

90°  plies.  Of  particular  interest  is  that  type  I  wet  and 

FPF 

stresses  predicted  at  are  nearly  equal.  A  picture  of  a  replica 
FPF 

taken  at  for  the  type  I  wet  laminate  is  shown  in  Figure  17  and 
demonstrates  that  both  transverse  cracking  and  delaminations  occur 
simultaneously  which  supports  the  approximation  of  equating  transverse 
strengths  in  the  lamina  2-3  plane  for  wet  90°  plies. 

In  summary,  the  variations  in  and  stresses  near  the  laminate 
edge  are  influenced  by  stacking  sequence  and  environmental  conditioning. 
Variations  of  and  stresses  in  the  90°  plies  near  the  edge  of  wet 
or  dry  type  I  and  type  II  laminates  cause  the  differences  in  experimental 
laminate  loads  required  to  initiate  cracking  in  the  90°  plies.  Although 
the  exact  magnitude  of  and  which  initiate  transverse  cracking  and 
delaminations  at  the  laminate  edge  are  not  known,  the  model  developed 
in  this  section  can  relate  the  experimental  laminate  FPF  loads  by 
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calculating  an  effective  transverse  edge  strength  which  is  assumed 
constant  within  a  small  region  near  the  laminate  edge.  Trends  and 
differences  of  the  experimental  laminate  FPF  loads  are  verified  by 
the  model  for  the  wet  and  dry  conditions  of  type  I  and  type  II  quas 
isotropic  laminates. 


V.  SUMMARY 


The  objective  of  this  Investigation  was  to  demonstrate  the  effect 
of  moisture,  residual  thermal  curing  stresses  and  mechanical  load  on 
the  damage  development  in  quasi-isotropic  laminates-  In  particular 
this  investigation  was  concerned  with  demonstrating  how  the  maximum 
moisture  absorbed  (wet)  in  type  I  and  type  II  laminates,  fabricated 
from  T300/5208  graphite/epoxy,  significantly  alters  the  dry  stress 
state  and  subsequent  damage  development  along  the  laminate  edge. 

Emphasis  is  placed  on  using  improved  values  for  v/et,  dry,  and  out- 
of-plane  elastic  properties  since  these  properties  are  required  to  pre¬ 
dict  the  damage  free  stress  state  at  the  laminate  edge.  Classical  lam¬ 
inate  theory  and  a  previously  developed  finite  eleii^ent  model  (FEM)  were 
used  to  predict  stress  states  prior  to  the  first  formation  of  damage. 
Crack  patterns  characteristic  of  the  laminate  in  a  wet  or  dry  condition 
were  also  predicted  using  a  previously  developed  shear  lag  model. 

Development  of  edge  damage  was  recorded  by  using  an  established 
replicating  technique  which  transfers  an  image  of  edge  damage  on  to  a 
thin  acetate  sheet.  Replicas  taken  during  the  test  can  be  immediately 
viewed  on  a  standard  microfiche  card  reader  which  allows  the 
experimentalist  to  interpret  the  edge  damage  and  interact  with  the 
test. 

The  effect  of  moisture  on  [Og],  [90g],  type  I  and  type  II  test 
specimens  is  suniriarized  in  Tables  14,  15,  and  16.  In  general  moisture 
tends  to  reduce  the  m  Julus,  strength,  and  scatter  of  strength  for 
all  specimens  tested.  The  difference  between  the  type  I  and  type  II 
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wet  strengths  is  also  observed  to  be  less  than  the  difference  between 
type  I  and  type  II  dry  strengths.  This  decrease  in  the  difference 
between  type  I  and  type  II  laminate  strengths  is  compared  as  shovm  in 
Figure  16  with  the  decrease  in  the  scatter  of  [Og]  strengths  when 
moisture  is  absorbed.  The  addition  of  moisture  in  the  [90g]  tests 
reduces  the  first  fracture  strength  from  a  dry  value  of  6.47  ksi  to 
a  wet  value  of  4.74  ksi  and  eliminates  the  dry  strength  equilibrium. 

The  absorbtion  of  moisture  also  causes  swelling  transverse  to 
the  fiber  direction  such  that  in  a  type  I  or  type  II  laminate  the 
residual  transverse  stress  in  the  90°  ply  changes  from  a  dry  value  of 
+3.43  ksi  to  a  wet  value  of  -1.96  ksi,  as  predicted  by  laminate  theory, 
when  1.2  percent  moisture  is  absorbed.  The  laminate  load  which  results 
in  the  first  90°  ply  failure  (FPF)  can  be  calculated  by  using  the 
maximum  stress  theory  which  assumes  FPF  occurs  when  the  laminate  theory 
prediction  of  transverse  stress  in  the  90°  ply  reaches  the  [90g]  test 
specimen  strength.  As  a  result  of  the  wet  residual  compressive  stress 
in  the  90°  ply,  the  predicted  laminate  load  required  to  cause  FPF  is 
increased  from  a  dry  value  of  789  Ib/in  to  a  wet  value  of  1900  Ib/in 
even  though  the  90°  ply  strength  is  reduced  due  to  moisture. 

The  first  occurrence  of  a  90°  ply  transverse  crack  is  a  single 
event  after  which  the  90°  ply  continues  to  carry  a  portion  of  the  FPF 
laminate  load.  As  the  laminate  load  is  increased  the  failure  process 
continues  due  to  the  formation  of  transverse  cracks  until  an  equilibrium 
spacing  betv/een  cracks  is  achieved.  When  equilibrium  spacing  is  ob¬ 
tained  the  90°  ply  no  longer  carries  the  stress  due  to  the  applied 
laminate  load. 


Moisture  was  also  shown  to  significantly  alter  the  dry  type  I 
laminate  edge  and  0^  stresses  predicted  by  the  FEH  such  that 
delamination  would  occur  at  a  lower  laminate  load  due  to  an  increase 
of  0^  when  moisture  is  absorbed.  Similarly  FPF  would  occur  at  a  higher 
laminate  load  due  to  a  predicted  decrease  in  0^  when  moisture  is 
absorbed.  Using  the  replica  technique  the  first  occurrence  of  wet  FPF 
was  observed  within  the  load  interval  (1100-1200)  Ib/in  which  is 
larger  than  the  FPF  laminate  load  of  (500-600)  Ib/in  which  was  observed 
for  the  dry  case.  Also  the  laminate  load  required  to  cause  delarnina- 
tions  was  observed  to  decrease  from  (1800-2000)  Ib/in  for  the  dry  case 
to  (1100-1200)  Ib/in  for  the  wet  case.  As  a  result  of  moisture 
absorbtion,  transverse  cracks  and  delaminations  were  observed  to  occur 
simultaneously  in  the  90°  plies  of  a  type  I  laminate  when  the  laminate 
load  reached  (1100-1200)  Ib/in. 

For  wet  and  dry  type  I  and  type  II  laminates,  FPF  laminate  loads 
were  observed  to  occur  at  lower  values  than  the  FPF  laminate  loads 
predicted  using  the  maximum  stress  theory  which  considers  only  the 
value  of  the  0^  stress  predicted  by  laminate  theory.  An  improved 
estimate  of  FPF  laminate  loads  should  include  the  and  edge 
stresses,  calculated  by  an  FEM  using  wet  and  dry  elastic  properties, 
and  interlaminar  strengths. 

For  wet  or  dry  laminates  there  is  a  difference  between  the  type 
I  and  type  II  FPF  laminate  loads.  The  difference  between  the  type  I 
and  type  II  wet  FPF  laminate  loads  was  larger  than  the  difference 
between  the  type  I  and  type  II  dry  FPF  laminate  loads.  A  model  was 


developed  v/hich  predicted  these  differences.  The  model  demonstrated 
that  differences  in  the  observed  FPF  laminate  loads  were  accounted  for 
by  the  predictable  changes  in  the  FEM  calculated  values  of  and 
edge  stresses  v/hich  were  the  result  of  changes  in  stacking  sequence 
together  with  wet  or  dry  environmental  conditioning. 

The  absorbtion  of  moisture  was  observed  and  predicted  to  have 
less  than  a  10  percent  increase  in  the  crack  spacing  between  the  90° 
ply  transverse  cracks  in  the  type  II  laminates.  For  the  type  I 
laminate  the  addition  of  moisture  nearly  doubles  the  distance  betv/een 
90°  ply  transverse  cracks  as  a  result  of  delaminations  occurrif'j 
simultaneously  and  interacting  with  the  transverse  cracks. 

Although  moisture  was  shov;n  to  significantly  alter  the  first 
formation  of  damage  in  the  90°  plies,  the  fully  developed  crack  patterns 
prior  to  fracture,  which  develop  from  static  or  cyclic  loads,  were  not 
significantly  altered  by  moisture.  Consequently,  the  difference 
between  the  redistributed  stresses  in  damaged  v/et  and  dry  laminates 
prior  to  fracture  will  be  small;  and  as  a  result,  these  differences 
will  have  a  negligable  effect  on  the  laminate  strength.  This  obser¬ 
vation  was  experimentally  verified  as  small  differences  between  wet 
and  dry  laminate  residual  or  static  strengths. 


VI.  CONCLUSIONS 


In  genera ‘i  moisture  tends  to  reduce  moduli,  strengths,  and  scatter 
of  strength  for  [Qg],  [90g],  [0/±45/90]^  and  [0/90/±45]^  specimens 
fabricated  from  T300/5208  graphite-epoxy.  Moisture  also  reduces  the 
difference  between  type  I,  [0/±45/90]^,  and  type  II,  [0/90/±45]^, 
laminate  strengths.  The  dry  equilibrium  strength,  which  is  experi¬ 
mentally  observed  for  the  dry  [90g]  tests,  is  eliminated  when  moisture 
is  absorbed. 

The  first  formation  of  edge  damage  such  as  first  90°  ply  failures 
(FPF)  and  delaminations  can  be  recorded  by  using  an  established  replica 
technique.  Damage  recorded  on  replicas  taken  during  the  test  can  be 
immediately  viewed  on  a  standard  microfiche  card  reader  such  that  the 
development  of  edge  damage  can  be  conveniently  interpreted  during  the 
test. 

Moisture  was  observed  to  significantly  alter  the  dry  type  I  laminate 
edge  residual  stress  state  predicted  by  a  FEM  such  that  delaminations 
would  occur  at  a  lower  laminate  load  a^'d  FPF  would  occur  at  a  higher 
laminate  load.  These  predicted  trends  due  to  moisture  absorbtion  were 
experimentally  verified  using  the  replica  technique. 

For  the  type  I  wet  laminate  delaminations  were  observed  to  occur 
simultaneously  and  interact  with  transverse  cracks  in  the  90°  plies 
such  that  the  equilibrium  spacing  v/as  twice  the  value  observed  fcr  the 
type  I  dry  laminate. 

The  absorbtion  of  moisture  was  observed  and  predicted  to  have  less 


than  a  10  percent  increase  in  the  equilibrium  crack  spacing  in  the  90° 
plies  of  the  type  II  laminate. 

By  using  the  model  developed  in  this  investigation  it  is  possible 
to  predict  changes  in  FPF  laminate  loads  due  to  differences  in  stacking 
sequence  together  with  v/et  or  dry  environmental  conditioning.  These 
differences  in  FPF  laminate  loads  were  accounted  for  by  the  predictable 
change  in  the  FEM  calculated  values  of  and  edge  stresses  due  to 
changes  in  stacking  sequence  together  with  wet  or  dry  environmental 
conditioning. 

Although  moisture  content  affects  the  load  at  which  damage  initi¬ 
ates  in  graphite  epoxy  laminates»  the  complete  damage  state  which 
develops  from  static  and  cyclic  loads  prior  to  fracture  is  a  charac¬ 
teristic  of  the  laminate  stacking  sequence  and  is  not  a  function  of 
loading  history  (monotonic  or  cyclic  loads)  and  environmental  condi¬ 
tioning  (wet  or  dry).  For  the  laminates  and  conditions  examined  in 
this  investigation,  the  experimental  data  show  that  the  tensile  strength 
of  monotonically  loaded  specimens  and  the  residual  tensile  strength  of 
cyclically  loaded,  fully  damaged,  specimens  are  dependent  on  stacking 
sequence  and  are  independent  of  the  hygro-mechanical  history  of  speci¬ 
mens  with  the  same  stacking  sequence.  The  results  suggest  that  strength 
of  a  composite  laminate  is  not  influenced  by  the  details  of  individual 
damage  events  but  rather  is  dependent  on  the  collective  form  of  the 
various  damage  details  as  described  by  the  concept  of  a  damage  state 
which  is  a  ’aminate  property  and  how  the  damage  state  affects  the 
strength  state  of  a  laminate. 
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APPENDIX  A 

TENSOR  POLYNOMIAL  FAILURE  CRITERION 


The  tensor  polynomial  failure  criterion  in  the  contracted  tensor 
notation  (for  an  orthotropic  material  in  the  principal  material 
directions)  has  the  form 


Fi  01+F2  1 

4.r  2,p  2  ,p  2  .p  2 

^33'’3  ^44'^23  ’^55‘^1 3  *^66^1 2 

+2Fi2Cia2+2F^30^C3+2F23C2a3  =  1 


(A.l) 


where  the  F.  and  F..  terms  are  as  previously  defined  in  Table  15. 

In  the  xyz  (laminate)  coordinate  system,  the  tensor  polynomial  failure 
criterion  transforms  (from  the  1-2  to  x-y  by  anticlockwise  rotation  of 
+0)  into 


F]0/F-ay+F'a^+F-axy+F'io2 

+''22'^y'^''33°z'^''44^yz'^''55^L 

^^^36^z\y*^^A5^yz\z^^^\2^x^y 

^2n3‘^x^z'^2f'23'^y"z  =  '' 


(A.2) 


where  the  F'  terms,  as  functions  of  the  unprimed  F's  and  e,  are  as 
follows  (m  =  cose,  n  =  sine) 


FI  =  ni^F,+n^F. 

I  1  C 

F^  =  n^F^+m^F^ 
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=  -2mn(F,-F2) 

F^l  =  m^F^^+m^n^(Fgg+2F^2)‘^"S2 

F^2  =  n'‘F^^+m^n^(Fgg+2F^2^^'"^''22 

"  *^33 

F^5  -  n2F4yF,3 

F^g  =  -mn[2(m^F^^-n^F22)-{ni^-f'^){2Fi2'^''66^^ 

F^g  =  -inn[2(n2F^^yF22)+{ni^-n^)(2F32+Fgg)] 

F'g  =  -mn{Fi3-F23) 

F45  =  "in(F44-Fgg) 

Ft2  “  fi^f'^{F]1+F22-Fgg)+(ni +n  )Fi2 
*"13  ""  "’^''l3'^"^^23 

Thpse  drs  transformations  from  th9  right  handed  1-2  coordinate  system 
into  another  right  hand  coordinate  system  obtained  by  an  anticlockwise 
rotation  of  0®  about  the  3  axis.  If  a  ply  is  oriented  at  +e®  from  the 
laminate  axis,  the  F..  are  obtained  by  using  the  above  equations  with 
the  sines  and  cosines  of  -0®. 
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APPENDIX  B 

ELASTIC  LAMINA  PROPERTIES 

Equations  developed  by  Hashin  [27]  are  listed  in  terns  of  fiber 
and  matrix  properties  denoted  by  subscripts  f  and  b,  respectively. 
The  fiber  properties  are  further  denoted  by  subscripts  L  and  T  which 
correspond  to  the  longitudinal  and  transverse  fiber  coordinates  de¬ 
fined  in  Figure  1.  Fiber  and  matrix  volume  fractions  are  denoted  by 
and  V^. 

The  upper  and  lower  bounds  for  the  plane  strain  bulk  modulus  in 
the  2-3  plane  shown  in  Figure  1  are  given  as 


f-l*  ’'f 

=  *<6"  Y--'  V,- 


(B.l) 


.{+)  -  K  +  ^ 

^23  ^f  V 

TT  1  ,  V 

K.-K,  K.  +G, 

^  ^T  ^TT  ^TT 


Where  *  indicates  the  equation  used  for  curve  fitting  and  the  (+)  and 
(-)  signs  indicate  upper  and  lower  bounds  respectively. 

The  remaining  equations  are  written  using  similar  notation. 


(a)  (B.2) 
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